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This  report  is  the  second  volume  of  a two-volume  final  report 
prepared  under  Contract  No.  F08606-74-C-0011  for  the  period  1 July  1975 
through  30  September  1976.  The  work  reported  on  involved  making  modi- 
fications to  the  computer  programs  previously  developed  for  use  in 
studying  the  atmospheric  transport  of  the  thermal  radiation  produced 
by  nuclear  sources  and  in  performing  parametric  studies  with  these 
programs. 


The  POLO  Monte  Carlo  procedure,  which  computes  the  transport  of 
thermal  radiation  from  point  sources  through  an  altitude-dependent 
atmosphere,  was  modified  to  utilize  the  gaseous  absorption  cross  section 
data  for  CO2,  H2O,  0^,  N2O,  CO,  CH^,  and  O2  that  had  been  compiled  by 
McClatchey  at  the  Air  Force  Geophysics  Laboratory.  The  AFGL  gaseous 
absorption  cross  section  data  are  available  on  magnetic  tape  for  use 
in  POLO.  The  gaseous  absorption  cross  section  data  have  a resolution 
of  Av  = 20  cm  ^ at  5 cm  ^ spacing  between  10  cm  ^ and  13155  cm  The 
modified  POLO  procedure  was  used  to  compute  direct  transmission  data 
at  a satellite  for  sources  positioned  within  two  model  atmospheres 
for  comparison  with  data  computed  with  the  L0VrTRAN3  code  for  the  same 
atmosphere.  The  good  comparison  obtained  between  the  results  of  the 
POLO  and  LOWTRAN  calculations  for  the  two  different  model  atmospheres 
verified  that  the  cilculational  method  used  in  POLO  properly  accounts 
for  gaseous  absorption.  In  addition,  POLO  can  be  used  to  compute  trans- 
missions less  than  0.001,  which  is  the  smallest  value  that  can  be 
computed  with  use  of  the  L0WTRAN3  procedure. 

The  POLO  procedure  was  used  to  study  the  transmission  of  infra- 
red radiation  from  point  isotropic  sources  to  satellite  positions  for 
wavelengths  between  2.6525  and  3.0030  pm.  The  gaseous  absorption  in 
that  wavelength  band  is  due  principally  to  H2O  and  CO2.  The  trans- 
mission calculations  were  made  with  the  POLO  procedure  for  both  a 
tropical  and  a mid-latitude  summer  atmosphere.  The  satellite  altitude 


was  t^en  to  be  35800  km  and  detectors  were  positioned  at  polar  position 
angles  (polar  angle  between  a radial  through  the  source  and  a radial 
through  the  satellite's  position)  of  0,  30,  60,  78,  and  82  degrees. 

The  source  altitudes  considered  were  0.1,  0.5,  2,  5,  and  10  km.  The 
point  isotropic  source  strength  was  taken  to  be  one  photon  uniformly 
distributed  across  the  source  wavelength  interval.  The  detector  filter 
function  was  taken  to  be  1.0  for  all  wavelengths  in  the  source  wave- 
length interval.  From  the  results  of  the  POLO  calculations,  it  was 
found  that  for  source  heights  above  3 km  the  direct  transmission  for 
the  mid-latitude  atmosphere  was  higher  than  that  computed  for  the  tropi- 
cal atmosphere.  For  source  heights  below  3 km  the  direct  transmission 
for  the  mid-latitude  atmosphere  was  less  than  that  computed  for  the 
tropical  atmosphere.  When  the  ground  albedo  was  0.0,  the  fraction  of 
the  total  transmittance  contributed  by  atmospheric  scattering  was  found 
to  Increase  with  an  increase  in  the  satellite  polar  position  angle  and 
with  a decrease  in  the  source  altitude.  When  the  ground  albedo  was  0.8, 
the  fraction  of  the  total  transmittance  contributed  by  atmospheric 
scattering  and  ground  reflection  decreases  with  an  Increase  in  the 
satellite  polar  position  angle  for  angles  up  to  about  65“  and  then 
increases  with  a further  Increase  in  the  satellite  polar  position  angle. 
The  POLO  calculations  show,  for  point  sources  emitting  in  the  2.7  ym 
wavelength  band,  that  the  air-scattered  and  ground-reflected  radiation 
can  be  an  important  component  of  the  radiation  field  measured  by 
satelllte-bome  sensors  in  that  wavelength  band. 

The  POLO  procedure  was  used  to  study  the  effect  of  source  aniso- 
tropy on  the  time-dependent  thermal  radiation  signal  at  satellite  posi- 
tions. The  calculations  were  run  for  a satellite  altitude  of  20187  km, 
source  altitudes  of  0.5,  3,  and  10  km  and  for  several  receiver  look 
angles  between  0“  and  90“ . The  sources  were  positioned  in  model  atmo- 
spheres having  sea-level  meteorological  ranges  of  3,  10,  and  25  km. 

The  time-dependent  transmission  data  generated  with  the  POLO  procedure 
was  convoluted  with  source  data  generated  with  the  RADFLO  code  for  a 
2-KT  event  and  detector  response  data  for  a wide-band  silicon  detector. 
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The  POLY  procedure  was  used  in  performing  the  convolution  calculations. 
The  times  at  which  the  power  is  one-half  and  one-tenth  of  that  at  first 
maximum  for  the  unattenuated  source  power  time  distribution  was  deter- 
mined to  be  51.72  microseconds  and  39.09  microseconds,  respectively. 

The  effect  of  atmospheric  scattering  and  ground  reflection  is  to  shift 
these  times  to  later  times.  Not  only  did  the  anisotropy  of  the  source 
affect  the  magnitude  of  the  optical  signal  at  the  detector  position, 
but  the  time  at  which  first  maximum  is  reached  was  also  shifted  to 
later  times  when  the  source  was  peaked  in  the  backward  direction. 

The  source  anisotrophy  has  a small  effect  on  the  times  of  one-half  and 
one-tenth  maximum  power.  For  small  look  angles,  the  peaking  of  the 
source  in  the  backward  direction  shifts  the  time  of  one-half  maximum 
power  by  approximately  0.85  microseconds  and  the  time  of  one-tenth 
maximum  power  by  approximately  0.52  microseconds.  For  the  90°  look 
angle,  the  average  shift  in  the  time  of  one-half  maximum  power  is  2.66 
microseconds  and  the  average  shift  in  the  time  of  one-tenth  first  maxi- 
mum power  is  1.13  microseconds.  However,  for  those  cases  where  the 
source  was  at  500  meters  altitude,  the  shift  in  time  of  one-half  maxi- 
mum power  is  as  much  as  7.48  microseconds  and  the  shift  in  the  time  of 
one-tenth  maximum  power  is  as  much  as  3.78  microseconds. 

The  POLO  procedure  was  used  to  extend  the  atmospheric  trans- 
mission data  base  previously  generated  for  meteorological  ranges  of 
3,  10,  and  25  km  to  a meteorolgoical  range  of  40  km.  The  time-dependent 
atmospheric  scattering  data  for  the  40-km  meteorological  range  atmo- 
sphere was  computed  for  wavelengths  of  0.4278,  0.5,  0.6,  0.75,  and 
1.07  vua,  source  altitudes  of  0.1,  0.5,  3,  and  10  km  and  receiver  look 
angles  of  0,  30,  45,  60,  70,  75,  80,  and  90  degrees.  Curve-fit  coeffi- 
cients for  use  in  the  TMTAU  program  were  generated  for  ground  albedo 
values  of  0.0,  0.2,  0.5,  and  0.8.  The  curve-fit  coefficient  data  are 
stored  on  Tape  No.  001706  at  AFTAC. 


The  POLO  data  for  the  3-,  10- , and  25-km  meteorological  range 
atmosphere  were  previously  computed  for  only  three  source  wavelengths. 
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0.4278,  0.75,  and  1.07  ym.  A procedure  denoted  as  POLICE  was  developed 
for  use  in  generating  curve-fit  coefficients  for  wavelengths  other  than 
the  three  wavelengths  used  in  the  POLO  calculations.  A modified  version 
of  TMTAU,  denoted  as  TMTRFOl  was  developed  to  use  curve-fit  coefficients 
for  more  than  three  source  wavelengths  in  evaluating  the  convolution 
integral  giving  the  time-dependent  power  at  receiver  positions.  In 
addition,  the  POLY  procedure  was  modified  to  accept  source  data  generated 
by  the  RADFLO  procedure  and  to  interpolate  the  atmospheric  transport  data 
over  wavelength  to  obtain  transport  data  at  the  wavelengths  used  in 
RADFLO  to  compute  the  wavelength  and  time  dependence  of  the  source 
emission.  The  TMTAU,  TMTRFOl,  and  POLY  procedures  were  checked  out 
against  an  analytic  solution  of  the  convolution  Integral  to  illustrate 
that  these  procedures  correctly  evaluate  the  convolution  Integral. 

A study  was  performed  with  the  POLO  procedure  to  determine  if 
the  time-dependent  atmospheric  transport  data  for  a given  satellite 
altitude  could  be  converted  to  give  data  for  other  satellite  altitudes. 
The  results  of  the  POLO  calculations  showed  that  direct  and  time- 
dependent  scattered  fluxes  for  satellite  altitudes  between  20186.8  km 
and  93621.626  km  could  be  obtained  from  data  previously  generated  with 
POLO  for  a satellite  altitude  of  35800  km  with  an  error  of  less  than 
1.0%.  It  was  found  that  the  correction  factors  for  converting  the  POLO 
data  for  a satellite  altitude  of  35800  km  to  satellite  altitudes 
between  20186.8  km  and  93621.626km  varied  significantly  with  the 
receiver  look  angle.  There  was  no  significant  variation  in  the  correc- 
tion factor  for  a given  look  angle  with  source  height. 

A number  of  POLO  problems  were  run  to  investigate  the  effect 
of  clouds  on  the  time-dependent  response  at  a satellite-borne  receiver. 
POLO  problems  were  run  for  10-km  meteorological  range  Haze  C atmosphere 
containing  either  a 20-mean-free-path-thick  altostratus  cloud  between 
3 and  4 km  altitude  or  a 14-mean-free-path- thick  stratus  cloud  at 
altitudes  between  0.5  and  1.0  km  altitude.  Point  isotropic  0.4278, 

0.75,  and  1.07  pm  wavelength  sources  were  positioned  at  altitudes  of 
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0.1,  0.5,  3,  and  10  km.  The  clouds  were  treated  as  scattering  media 
and  the  ground  surface  was  considered  to  be  a Lambert  reflector.  The 
satellite  altltvide  was  taken  to  be  35800  km.  The  POLO-generated  data 
were  convoluted  with  a RADFLO-generated  source  term  for  a 2-KT  event 
to  investigate  the  effect  of  the  cloud  layer  on  the  times  at  which 
first  maximum  power  and  one-tenth  first  maximum  power  occurs.  When 
the  source  is  positioned  under  a cloud  layer,  the  retarded  time  of 
first  maximum  power  tends  to  Increase  with  look  angle.  The  retarded 
time  at  which  one-tenth  of  first  maximum  power  occurs  also  tends  to 
Increase  with  look  angle.  When  the  source  is  above  either  cloud  layer 
at  10-km  altitude,  the  retarded  time  at  which  first  maximum  occurs 
tends  to  decrease  with  an  increase  in  the  look  angle  and  the  time  at 
which  one-tenth  of  first  maximum  power  occurs  remains  nearly  a constant 
for  all  look  angles. 

Attempts  to  use  POLO  to  compute  the  time-  and  angle-dependent 

distributions  of  the  light  escaping  from  the  top  of  a cloud  located 

above  a point  isotropic  source  resulted  in  data  with  large  statistical 

errors.  These  large  statistical  errors  are  believed  to  be  caused 

-2 

by  the  use  of  a point-estimating  function  that  has  an  R dependence 
on  the  source-receiver  distance  R.  The  large  statistical  variation 
in  the  POLO  data  tend  to  become  insignificant  when  the  receivers 
are  positioned  at  satellite  altitude.  An  effort  has  been  initialed 
to  develop  a version  of  POLO  that  uses  a surface-receiver  estimating 
function  Instead  of  a point-receiver  estimating  function  to  deter- 
mine the  radiation  intensity  at  receiver  location  at  the  top  of  the 
cloud. 

A study  was  initiated  to  Investigate  the  possible  application 
of  discrete  ordinates  methods  of  radiation  transport  to  study  time- 
dependent  radiation  transport  through  clouds.  The  use  of  discrete 
ordinates  methods  requires  that  the  phase  function  be  expanded  in 
terms  of  Legendre  polynomials.  A method  for  computing  the  coefficients 


vil 


for  the  Legendre  polynomial  expansion  from  the  use  of  Mie  theory  was 
developed  but  it  was  found  that  the  number  of  coefficients  required  to 
represent  the  phase  function  is  on  the  order  of  2a  + 2 where  a = 2iTr/X 
is  the  aerosol  size  parameter.  Clouds  typically  have  aerosols  up  to 
20  to  30  microns  in  radius.  For  0.5  ym  wavelength  light,  the  size 
parameter  could  be  as  large  as  375,  thus  requiring  a Legendre  poly- 
nomial of  the  order  of  752. 
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I . INTRODUCTION 


This  report  is  the  second  volume  of  a two-volume  Final  Report 
prepared  under  Contract  No.  F08606-74-C-0011  for  the  period  1 July  1975 
through  30  September  1976.  The  work  reported  on  involved  making  modi- 
fications to  the  computer  programs  previously  developed  for  use  in 
studying  the  atmospheric  transport  of  the  thermal  radiation  produced 
by  nuclear  sources.  The  work  also  involved  applying  these  computer 
procedures  to  a number  of  different  thermal  radiation  transport  problems. 

The  POLO  Monte  Carlo  procedure  (Ref.  1),  which  computes  radiation 
transport  for  point  sources  in  a spherical- shell  atmosphere,  was  modified 
in  order  to  utilize  the  gaseous  absorption  cross-section  data  compiled 
by  McClatchey  (Ref.  2)  at  AFGL.  These  new  cross-section  data  allow  the 
POLO  procedure  to  treat  absorption  by  a larger  number  of  atmospheric 
gases  and  to  treat  gaseous  absorption  within  finite-sized  wavelength 
intervals  in  the  infrared  wavelengths  Instead  of  for  discrete  wavelengths. 
The  POLO  procedure  was  also  modified  so  that  it  could  compute  and  print 
out  the  intensity  of  the  scattered  radiation  at  point  receivers  as  a 
function  of  a polar  and  an  azimuthal  angle  as  well  as  time.  The  modi- 
fications made  to  POLO  for  the  purpose  of  incorporating  the  AFGL  gaseous 
absorption  cross-section  data  and  printing  out  the  polar  and  azimuthal 
angle  distributions  of  the  scattered  Intensities  as  a function  of  time 
are  described  in  Section  II. 

The  modified  version  of  POLO  was  used  to  do  a parametric  study 
of  the  effects  of  atmospheric  scattering  on  the  transmission  of  the 
radiation  emitted  by  point  isotropic  sources  in  the  2.7  ym  band.  The 
atmospheric  transmission  to  a satellite-based  receiver  position  was 
evaluated  as  a function  of  the  satellite  look  angle  for  a satellite 
at  35,800  km  altitude  and  point  isotropic  sources  at  altitudes  of  0.1, 

0.5,  2,  5,  and  10  km  in  both  tropical  and  midlatitude  summer  model 
atmospheres.  The  atmospheric  models  and  the  results  obtained  from  the 
POLO  calculations  are  discussed  in  Section  III. 
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A study  was  conducted,  using  the  POLO  and  POLY  (Ref.  1)  pro- 
cedures, to  determine  the  effect  of  source  anisotropy  on  the  time  dis- 
tribution of  the  thermal  radiation  signal  that  would  be  seen  by  a wide- 
band silicon  receiver  positioned  on  a satellite  at  20,l87-km  altitude. 

A description  of  the  results  obtained  for  two  different  anisotropic 
sources  and  for  a point  isotropic  source  is  given  in  Section  IV. 

Section  V describes  a set  of  calculated  atmospheric  transport 
data  that  were  generated  with  the  POLO  procedure  for  a 40-km  meteoro- 
logical range  atmosphere.  The  receivers  were  at  synchronous  satellite 
altitude.  The  time-dependent  air-scattering  data  generated  with  POLO 
were  curve-fitted  for  use  with  the  TMTAU  procedure  (Ref.  3).  The 
curve-fit  data  are  also  discussed  in  Section  V. 

Section  VI  describes  the  calculatlonal  methods  developed  for 
convoluting  time-  and  spatial-dependent  source  and  atmospheric  trans- 
mission data.  Modifications  that  were  made  to  the  POLY  and  TMTAU  pro- 
cedures so  that  they  could  interpolate  the  POLO- genera ted  atmospheric 
data  or  the  curve-fit  coefficients  for  the  POLO  data  over  wavelength 
are  also  described.  The  accuracy  of  the  calculatlonal  methods  used  in 
the  modified  convolution  procedures  was  checked  through  comparison  with 
the  results  of  an  analytic  solution.  The  results  obtained  are  also 
discussed  in  Section  VI. 

The  effects  of  clouds  on  the  time-dependent  transmission  of 
the  thermal  radiation  emitted  by  nuclear  sources  was  studies  with  use 
of  the  POLO  and  POLY  procedure.  POLO  calculations  were  performed  for 
point  isotropic  sources  at  wavelengths  of  0.4278,  0.75,  and  1.07  pm,  and 
source  heights  of  0.1,  0.5,  3,  and  10  km.  The  detectors  were  positioned 
at  35,800-km  altitude  for  look  angles  of  0,  45,  75,  and  90  degrees.  Two 
model  atmospheres  were  used  in  the  POLO  calculations.  The  first  model 
atmosphere  used  was  a 10-km  meteorological  range  Haze  C atmosphere  with  a 
14  mean-free-path  thick  stratus  cloud  positioned  between  0.5-  and  1.0-km 
altitude.  The  second  atmosphere  used  was  a 10-km  meteorological  range 
Haze  C atmosphere  with  a 20  mean-free-path  thick  altostratus  cloud  posl- 


2 


tloned  between  3-  and  4-km  altitude.  The  POLO-computed  transmission 
date  were  convoluted  in  the  POLY  procedure  with  source  data  from  RADFLO 
for  a 2-KT  event.  The  results  of  both  the  POLO  calculations  of  the 
atmospheric  transmittance  and  the  POLY  calculations  of  the  thermal  power 
at  the  detector  for  retarded  times  up  to  140  psec  are  described  in 
Section  VII. 

The  POLO  procedure  was  used  to  compute  the  angular  and  time 
distributions  of  the  scattered  intensity  at  the  top  of  a cloud.  An 
analysis  of  those  calculations  showed  that  large  statistical  deviations 
occurred  in  the  calculational  data.  A discussion  of  these  results  and 
possible  reasons  for  the  large  statistical  errors  for  detectors  at 
cloud  top  is  also  given  in  Section  VII. 

There  have  been  several  occasions  where  calculated  air  trans- 
mission data  for  satellite  altitudes  other  than  that  of  35800  km  have 
been  needed.  There  currently  exists  a large  body  of  POLO-calculated 
atmospheric  transport  data  for  a satellite  altitude  of  35800  km.  A 
study  was  performed  utilizing  the  POLO  procedures  to  determine  the  time- 
dependent  scattered  fluxes  reaching  satellites  at  altitudes  of  20186.8, 
93261.626,  and  35800  km.  Ratios  of  the  calculated  scattered  fluxes 
at  altitudes  of  20186.8  and  93621.626  km  to  the  calculated  scattered 
flux  at  35800-km  altitude  were  found  to  agree  within  a few  percent 
with  the  similar  ratios  computed  for  the  direct  fluxes.  Because  of 
this  good  agreement,  it  is  suggested  that  the  method  described  in 
Section  VI  be  used  to  convert  scattered  fluxes  and  direct  fluxes  com- 
puted by  use  of  POLO  for  an  altitude  of  35800  km  to  any  other  satellite 
altitude. 

The  difficulties  in  computing  accurate  transmission  data  for 
detectors  positioned  at  cloud  top  with  the  POLO  program  resulted  in 
a study  on  the  possible  application  of  discrete  ordinates  techniques 
to  cloud  transmission  problems.  The  TDA  code  (Ref.  4),  a one- dimensional 
discrete  ordinates  code  was  obtained  from  the  Radiation  Shielding  Infor- 
mation Center.  It  soon  became  apparent  that  in  order  to  be  able  to 
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use  TDA,  or  any  other  time- dependent  discrete-ordinates  code,  to  study 
optical  radiation  transport  it  would  be  necessary  to  have  a calcula- 
tional  method  for  expanding  the  aerosol  phase  function  in  terms  of  a 
Legendre  polynomial  of  order  n.  A method  was  developed  to  accomplish 
this  expansion  from  data  computed  with  the  MIE-2  procedure.  The  cal- 
culatlonal  method  developed  to  compute  the  Legendre  coefficients  for 
Mie  scattering  is  discussed  in  Section  VIII. 
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II.  MODIFICATIONS  TO  THE  POLO  PROCEDURE 


During  the  previous  year,  several  modifications  were  made  to 
the  POLO  program  to  improve  the  application  of  the  program  to  problems 
of  Interest  to  AFTAC.  Currently,  three  versions  of  the  POLO  program 
are  operational  on  the  computer  at  AFTAC.  The  original  version  of  POLO, 
designated  TPOLOW,  is  still  operational  and  is  being  used  in  time- 
dependent  light-transport  studies  for  wavelengths  in  the  visible  range. 
This  version  of  the  POLO  program  produces  an  output  tape  in  the  format 
required  by  the  TMTAU  and  POLY  programs. 

In  a version  of  POLO,  designated  as  TPOLOG,  a new  method  for 
treating  gaseous  absorption  was  inc^/rporated  into  the  POLO  program.  The 
gaseous  absorption  data  compiled  by  McClatchey  (Ref.  2)  for  CO^,  H^O, 

O3,  N^O,  CO,  CH^,  and  0^  are  contained  on  a library  tape  which  is  ana- 
lyzed by  the  TPOLOG  program  to  select  the  cross-section  data  for  the 
desired  constituents  and  wavelength  intervals.  The  data  are  available 
on  tape  with  a resolution  of  Av  = 20  cm  ^ at  5 cm  ^ spacing  between 
10  cm  ^ and  13155  cm  (0.76-1000  Um)  . Therefore,  TPOLOG  treats 
gaseous  absorption  in  the  IR  for  a wavelength  band  rather  than  for 
a discrete  wavelength.  The  gaseous  absorption  transmission  model 
used  in  TPOLOG  with  the  McClatchey  data  is  described  in  Ref.  5. 

The  TPOLOG  program  was  further  modified  to  create  a version  of 
POLO  designated  as  TPOLOL.  TPOLOL  allows  for  a detector  with  polar 
angle,  azimuthal  angle,  and  time  resolution.  The  polar  angle  at  the 
detector  is  measured,  depending  on  the  option  selected,  either  from  the 
earth  radial  through  the  receiver  position  or  from  the  source-receiver 
axis.  An  effort  was  made  to  minimize  the  differences  in  the  utiliza- 
tion of  the  TPOLOG  and  TPOLOL  procedures.  TPOLOL  currently  exists  at 
AFTAC  on  disc  pack  TF2PAK  under  RRACODES  in  all  sections  as  TPOLOL. 

2.1  Comparison  of  TPOLOG  Calculations  with  L0WTRAN3  Data 

(U)  TPOLOG  calculations  of  the  direct  transmittance  in  the 
2,65  pm  to  3.0  pm  wavelength  range  for  a point  source  at  a height  of 


A 
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0.5  km  and  a detector  positioned  at  a satellite  altitude  of  35800  km 
were  run  for  comparison  with  L0WTRAN3  (Ref.  6)  calculations.  The  cal- 
culations were  made  for  two  model  atmospheres.  The  atmosphere  denoted 
as  Model  A had  less  water  vapor  content  than  did  the  atmosphere  denoted 
as  Model  B.  The  LOWTRAN3  code  was  run  with  the  same  atmospheric  models 
as  were  used  in  the  TPOLOG  calculations.  A comparison  of  the  direct- 
transmittance  data  obtained  from  both  calculations  for  the  two  atmo- 
spheric models  when  the  satellite  is  positioned  directly  overhead  is 
shown  in  Fig.  1 for  the  range  2.8249  ym  - 2.9940  ym  (3,540  cm  ^ and 
3,340  cm  respectively).  In  the  figure,  the  direct-transmittance 
data  obtained  from  L0WTRAN3  are  connected  with  straight  lines;  whereas 
the  TPOLOG-calculated  direct-transmittance  data  are  plotted  using  circles 
and  triangles.  The  differences  that  occur  for  the  individual  data  is 
believed  to  be  due  to  the  different  method  of  treating  absorption  in 
TPOLOG  and  in  L0WTRAN3.  TPOLOG  uses  an  interpolation  between  absorption 
cross  sections  as  a function  of  temperature  in  order  to  give  the  correct 
cross  section  at  the  specific  temperature  of  each  atmospheric  layer. 

The  transmission  formula  in  TPOLOG  (Ref.  5)  takes  into  account  the 
pressure  dependency.  The  absorption  coefficients  in  L0WTRAN3  are  input 
for  296“K  temperature  and  1 atm  pressure  and  are  then  adjusted  to  the 
correct  temperature  and  pressure  of  each  atmospheric  layer  (the  adjust- 
ment takes  into  account  the  change  of  the  line  widths  with  temperature, 
but  not  the  change  of  the  line  strengths) . 

Unfortunately,  the  transmission  functions  Incorporated  in 
L0WTRAN3  do  not  allow  the  use  of  the  program  to  calculate  transmissions 
that  are  less  than  0.001  for  most  of  the  species  considered.  L0WTRAN3 
transmission  data  in  the  range  2.65  ym  were  all  calculated  to  be  0.0 
and  therefore,  cannot  be  compared  with  the  TPOLOG  data  in  the  same 
wavelength  range. 

A comparison  of  the  direct  transmit tances  for  the  Model  B atmo- 
sphere that  were  calculated  by  both  programs  is  shown  in  Fig.  2 for 
satellite  polar  position  angles,  0^,  of  0°,  lO” , 60°,  and  78°.  The 
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comparison  is  good  for  transmittances  greater  than  0.001  and  for  all 

of  the  polar  position  angles,  with  the  exception  that  the  L0WTRAN3 

data  seem  to  be  slightly  lower  than  the  TPOLOG  data  for  0 = 60“  and 

o 

78“ . This  is  explained  by  the  fact  that  L0WTRAN3  considers  the  bending 
of  light  by  the  atmospheric  refraction;  whereas  the  version  of  TPOLOG 
that  was  used  to  calculate  the  infrared  data  does  not.  The  atmospheric 
refraction  tends  to  bend  the  light  rays  to  lower  atmospheric  layers, 
resulting  in  larger  pathlengths  and  lower  transmissions.  Since  the 
overall  comparison  between  the  L0WTRAN3  and  TPOLOG  data  is  good,  it  was, 
therefore,  decided  not  to  modify  L0WTRAN3  to  neglect  the  effect  of  the 
refractive  index  which  would  require  major  changes  in  L0WTRAN3. 

The  good  comparison  between  transmission  data  computed  with 
LOWTRAfLl  and  TPOLOG  verifies  the  method  used  in  TPOLOG  to  account  for 
gaseous  absorption.  The  method  used  in  TPOLOG  seems  to  have  the  advan- 
tage that  transmission  data  can  be  obtained  for  values  smaller  than 
0.001;  however,  no  information  is  available  on  the  accuracy  of  the 
computed  data  for  very  small  transmissions. 


III.  TPOLOG  CALCULATIONS  IN  THE  IR  FOR  SATELLITE  DETECTORS 
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TPOLOG  multiple-scattering  calculations  In  the  IR  for  detector 
systems  at  satellite  positions  have  been  performed  In  the  range  2.65  - 

3.00  ym.  This  wavelength  Is  of  particular  Interest  because  of  the  H2O 
and  CO2  absorption  band  In  that  region  of  the  spectrum.  Because  of 
the  high,  but  different,  concentration  of  these  species  at  low  alti- 
tudes, H2O  absorption  Is  dominant  In  the  atmospheric  layers  close  to 
the  ground  surface;  CO^  absorption  Is  dominant  above  5-kllometer  alti- 
tudes and  can  still  be  observed  at  altitudes  about  40  km  (Ref.  7). 

3.1  Model  Atmospheres 

Atmospheric  models  for  a tropical  and  a mid-latitude  summer 
atmosphere  were  chosen.  The  dependency  of  the  atmospheric  pressure  on 
altitude  for  these  models  was  taken  out  of  the  U.  S.  1966  Supplement 
Atmospheres  (Ref.  8).  Data  for  temperature  and  dew  point  below  the 
tropopause  were  supplied  by  the  contract  monitor.  The  temperature 
data  obtained  were  smoothed  and  then  extended  to  higher  altitudes 
using  the  tropical  and  mid-latitude  models  of  Ref.  8. 

Although  the  TPOLOG  program  allows  for  gaseous  absorption  for  a 

total  of  seven  species,  only  the  absorption  by  H^O,  CO^,  and  N2O  have 

to  be  considered  In  the  wavelength  range  of  Interest  (Ref.  9).  CO2 

and  NO  was  assumed  to  be  uniformly  mixed  In  the  entire  atmosphere 
^ -4  -7 

with  mixing  ratios  of  3X10  and  5X10  (Ref.  10),  respectively.  The 
H2O  densities  below  lO-km  altitude  were  calculated  from  the  measured 
dew  point  and  temperature,  using  the  USAF  SKEW  T,  log  p Diagram 
supplied  by  the  contract  monitor,  and  tables  given  In  Ref.  11.  The 
data  were  then  smoothed  and  combined  with  data  given  by  Selby  et  al. 
(Ref.  9)  for  both  the  tropical  and  mid-latitude  summer  model. 

The  dependency  of  the  aerosol  concentration  upon  altitude  was 
assumed  to  be  that  of  the  1970  measurements  taken  by  Elterman  (Ref.  12). 
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Aerosol  models  for  tropospheric  and  stratospheric  particles  were  taken 
from  Ref.  13  to  describe  the  aerosol  size  distribution  and  the  wave- 
length dependency  of  the  aerosol  scattering  coefficient  (see  also 
Ref.  7).  The  ground-level  meteorological  ranges  for  both  the  tropical 
and  mid-latitude  atmospheres  were  taken  to  be  25  km. 

3.2  Calculational  Results 

TPOLOG  calculations  for  the  2.7-mlcron  wavelength  band  were  run 
for  the  tropical  and  mid-latitude  atmospheric  models.  The  gaseous 
absorption  cross-section  data  compiled  by  McClatchey  at  AFCRL  (Ref.  2) 
were  used  in  the  TPOLOG  calculations.  The  satellite  altitude  was  taken 
to  be  35,800  km  and  detectors  were  positioned  at  polar  position  angles 
(polar  angle  between  a radial  through  the  source  and  a radial  through 
the  satellite's  position)  of  0,  30,  60,  78,  and  82  degrees.  The  source 
altitudes  considered  were  0.1,  0.5,  2,  5,  and  10  km. 

The  TPOLOG  calculations  were  run  for  23  wavelengths  spaced 
between  2.6525  pm  and  3.0030  um.  The  scattering  and  continuum  absorp- 
tion cross  sections  were  input  for  three  wavelength  intervals;  the 
first  and  second  intervals  each  contained  seven  wavelengths  and  the 
third  interval  contained  nine  wavelengths.  For  the  2.7-ym  wavelength 
band,  the  atmospheric  gases  that  contribute  to  the  molecular  absorption 
were  water  vapor,  CO^,  and  N^O.  The  atmospheric  optical  thickness  for 
each  of  the  model  atmospheres,  excluding  gaseous  absorption,  were  as 
follows: 


WAVELENGTH 
INTERVAL  (um) 

OPTICAL  THICKNESS  EXCLUDING 

MOLECULAR  ABSORPTION 

MID-IATITUDE  ATMOSPHERE 

tropichl  atmosphere 

2.6525-2.7491 

4.033676X10"^ 

4.065740Xlu"^ 

2.7491-2.8490 

3.595484X10"^ 

3.625945X10"^ 

2.8490-3.0030 

3.334809X10"^ 

3.363596X10"^ 
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A significant  fraction  (45%  to  65%,  depending  upon  altitude)  of  the 
above  optical  thicknesses  resulted  from  aerosol  particle  absorption. 

The  TPOLOG  printed  output  gave  the  scattered  and  direct  flux 

at  each  receiver  position  as  a function  of  the  ground  albedo  for  each 

of  the  23  source  wavelengths.  The  ground  surface  was  assumed  to  be 

a Lambert  reflector.  Calculations  were  run  for  source  heights  of  0.1, 

0.5,  2,  5,  and  10  km.  The  receivers  were  positioned  at  an  altitude 

above  the  ground  surface  of  35,800  km  and  at  polar  position  angles 

0 of  0°,  30“,  60°,  78°,  and  82°. 
o 

The  calculated  direct  and  scattered  fluxes  were  folded  with 

a filter  function  and  a source  term  that  assumed  a uniform  emission 

with  wavelength  within  the  interval  being  considered.  The  results  were 

2 

then  multiplied  by  the  quantity  4ttR  where  R is  the  source-receiver  dis- 
tance to  obtain  the  direct  and  scattered  transmittances  to  the  satellite. 
The  total  source  emission  over  the  wavelength  interval  from  2.6525  pm 
to  3.0030  pm  was  assumed  to  be  1.0  photon.  The  results  shown  in  this 
report  are  for  a filter  function  that  had  a value  of  1.0  for  all 
wavelengths. 

The  results  obtained  from  the  TPOLOG  calculations  as  a function 
of  source  height  and  satellite  polar  position  angle  for  the  direct 
transmission  and  for  the  scattered  transmission  when  the  ground  albedo 
has  values  of  0.0  and  0.8  are  tabulated  in  Tables  I and  II  for  the  tropi- 
cal and  mid- latitude  atmospheres,  respectively.  From  the  data  shown  in 
these  tables,  it  is  seen  that  the  direct  transmittance  for  the  mid- 
latitude  atmosphere  is  higher  than  that  computed  for  the  tropical  atmo- 
sphere at  source  heights  less  than  about  3 km.  For  source  heights  above 
3 km,  the  direct  transmittances  for  the  mid- latitude  atmosphere  are  less 
than  those  computed  for  the  tropical  atmosphere.  Note  also  that  for  a 
satellite  polar  position  angle  of  82°  the  direct  transmittance  is  negli- 
gible compared  to  the  scattered  transmittance  when  the  source  height  is 
less  than  2 km. 

The  ratio  of  the  calculated  scattered  transmittance  to  the 
total  transmittance  is  shown  in  Fig.  3 for  the  tropical  atmosphere 


12 


TABLE  I.  TRANSMITTANCE  FOR  TROPICAL  ATMOSPHERE  VS  SOURCE  HEIGHT  AND 


Read  .36274-01  as  .36274X10 


TABLE  II.  TRANSMITTANCE  FOR  MID-LATITUDE  ATMOSPHERE  VS  SOURCE  HEIGHT 
AND  SATELLITE  POLAR  POSITION  ANGLE 


Read  .54960-01  as  .54960X10 


0 30  60  90 

SATELLITE  POLAR  POSITION  ANCLE  (deg) 


Fig.  3.  Ratio  of  the  Scattered-to- Total  Transmittance 
In  the  Tropical  Atmosphere  vs  Satellite  Polar 
Position  Angle  for  Source  Height  of  0.1,  0.5, 
2,  5,  and  10  km  and  a Ground  Albedo  of  0.0 
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as  a function  of  the  satellite  polar  position  angle  and  the  source 
height  for  a ground  albedo  value  of  0.0.  Similar  data  for  a ground 
albedo  value  of  0.8  are  shovm  in  Fig.  4.  Vfhen  the  ground  albedo  is 
0.0,  the  fraction  of  the  total  transmittance  contributed  by  atmospheric 
scattering  increases  with  an  increase  in  the  satellite  polar  position 
angle  and  with  a decrease  in  source  altitude.  When  the  ground  albedo 
is  0.8,  the  fraction  of  the  total  transmittance  contributed  by  atmo- 
spheric scattering  and  ground  reflection  decreases  with  an  increase 
in  the  satellite  polar  position  angle  up  to  about  65°  and  then  increases 
with  a further  increase  in  the  polar  position  angle.  As  seen  in  Figs. 

3 and  4,  the  importance  of  atmospheric  scattering  decreases  with  an 
increase  in  the  source  altitude.  It  should  be  noted  that  the  magnitude 
of  the  ground  albedo  affects  significantly  the  fraction  of  the  total 
radiation  at  the  satellite  that  is  contributed  by  the  atmospheric- 
scattered  and  ground-reflected  radiation.  The  TPOLOG  calculations  show 
that  the  magnitude  of  the  scattered  radiation  at  the  satellite  increases 
approximately  exponentially  with  an  increase  in  the  ground  albedo. 

Figure  5 shows  the  variation  of  the  ratio  of  the  scattered- 
to-total  transmittance  at  the  satellite  with  source  height  for  the  tropi- 
cal atmosphere.  Data  are  shown  in  Fig.  5 for  satellite  polar  position 
angles  of  0°  and  78°  and  for  ground  albedo  values  of  0.0  and  0.8.  For 
a satellite  polar  position  angle  of  0°,  the  magnitude  of  the  ratio  of 
the  scattered-to-total  transmittance  for  a given  source  height  is  highly 
dependent  on  the  magnitude  of  the  ground  albedo.  For  a polar  position 
angle  of  78°  and  for  a given  source  altitude,  the  ratio  of  the  scattered- 
to-total  transmittance  is  not  very  dependent  on  the  magnitude  of  the 
ground  albedo. 

The  TPOLOG  calculations  show,  for  point  sources  emitting  in 
the  2.7  ijm  wavelength  band,  that  the  air-scattered  and  ground-reflected 
radiation  can  be  an  important  component  of  the  data  measured  by  infrared 
sensors  on  satellites,  and  neglecting  that  component  in  an  analysis  of 
the  measured  data  can  lead  to  significant  error  in  estimating  the 
magnitude  of  the  source  that  produced  the  measured  data. 
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SCATTERED  TRANSMITTANCE /TOTAL  TRANSMITTANCE 


0 30  60  90 

SATELLITE  POLAR  POSITION  ANGLE  (deg) 


Fig.  4.  Ratio  of  the  Scattered-to-Total  Transmittance 
in  the  Tropical  Atmosphere  vs  Satellite  Polar 
Position  Angle  for  Source  Heights  of  0.1,  0.5, 
2,  5,  and  10  km  and  a Ground  Albedo  of  0.8 
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SCATTERED  TRANSMITTANCE /TOTAL  TRANSMITTANCE 


Fig.  5.  Ratio  of  the  Scattered-to-Total  Transmittance  in  the 
Tropical  Atmosphere  vs  Height  of  Source  for  Satellite 
Polar  Position  Angles  of  0“  and  78°  and  Ground  Albedo 
Values  of  C.O  and  0.8 
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IV.  POLO  CALCULATIONS  FOR  ANISOTROPIC  SOURCES 


A study  was  conducted  to  determine  the  effect  of  source  aniso- 
tropy on  the  time  distribution  of  the  optical  signal  received  by  a wide- 
band silicon  receiver  located  on  a satellite  positioned  20,187  kilometers 
above  the  earth's  surface.  Point  anisotropic  sources  were  considered 
which  emit  different  fractions  of  the  total  radiance  uniformly  in  the 
forward  and  backward  hemispheres.  The  forward  and  backward  hermlspheres 
are  separated  by  a vertical  plane  containing  the  source  point  that  is  per- 
pendicular to  the  plane  containing  the  receiver  locations.  The  forward 
hemisphere  is  the  hemisphere  on  the  receiver  side  of  the  plane.  The 
source-receiver  geometry  for  the  calculations  is  shown  in  Fig.  6.  Source 
altitudes  of  500  meters,  3 kilometers,  and  10  kilometers  were  considered 
in  model  atmospheres  having  meteorological  ranges  of  3,  10,  and  25  kilo- 
meters. Receivers  were  located  at  look  angles,  0^^,  of  0,  30,  45,  60,  70, 
75,  80,  and  90  degrees  for  the  sources  at  500  meters  altitude  and  at  0, 

30,  60,  70,  80,  and  90  degrees  for  the  sources  at  3 and  10  kilometers 
altitude. 

Atmospheric  transmission  calculations  for  each  source-receiver 
combination  in  each  of  the  three  model  atmospheres  were  performed  with 
the  POLO  program  for  source  wavelengths  of  0.4278,  0.75,  and  1.07  microns. 
The  direct  fluxes  that  were  obtained  for  the  three  different  0.75-micron 
wavelength  sources  located  at  500-meters  altitude  in  the  25-kilometer 
meteorological  range  atmosphere  are  compared  in  the  table  below.  The 
anlsotropic-1  source  emits  2/3  of  Its  radiation  in  the  forward  hemisphere 
and  1/3  in  the  backward  hemisphere.  The  anisotroplc-2  source  emits  2/3 
of  its  radiation  in  the  backward  hemisphere  and  1/3  in  the  forward  hemi- 
sphere. 

The  direct  flux  at  the  satellite  for  each  look  angle  and  source 
model  when  the  source  wavelength  is  0.75  um  is: 
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Fig.  6.  Source-Receiver  Geometry 
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DIRECT  FLUX  (photons  km  /source  photon) 


Cdeg) 

ISOTROPIC  SOURCE 

ANISOTROPIC-1 

ANISOTROPIC-2 

0 

.16346-9 

.16346-9 

.16346-9 

30 

.14801-9 

.19735-9 

.98673-10 

60 

.10504-9 

.14005-9 

.70027-10 

90 

.18777-14 

.24236-14 

.12118-14 

These  data  are  for  a source  height  of  0.5  km  In  the  25-km  meteorological 
range  atmosphere.  For  the  two  anisotropic  source  models  studies,  the 
direct  flux  at  6j^  = 0“  is  identical  to  that  for  an  Isotropic  source. 

A plot  of  the  variation  of  the  total  flux  at  the  detector  with 
a look  angle  for  a 0.75-M  wavelength  source  is  given  in  Fig.  7 for  a 
ground  albedo  of  0.0  and  in  Fig.  8 for  a ground  albedo  of  0.8.  The 
fluxes  in  Figs.  7 and  8 are  for  a source  altitude  of  0.5  km  and  a 25-km 
meteorological  range  atmosphere.  It  appears  that  the  possibility  of 
determining  whether  a source  is  anisotropic  or  not  with  the  use  of 
measured  data  from  at  least  two  satellites  is  the  greatest  for  look  angles 
between  40°  and  70°  and  when  the  surface  albedo  is  small.  The  two  satel- 
lites would,  of  course,  have  to  be  on  opposite  sides  of  the  source. 


The  atmospheric  transmission  computed  with  POLO  for  the  three 
wavelengths  were  used  in  the  POLY  procedure  to  predict  the  atmospheric 
transmission  for  14  wavelength  intervals  covering  the  range  from  0.4  to 
1.10  microns.  The  atmospheric  transmissions  for  each  of  the  14  wave- 
length intervals  were  convoluted  in  POLY  with  a time-dependent  source 
term  generated  with  the  RADFLO  program  for  a 2-KI  event.  Then  these 
data  were  weighted  in  POLY  with  a response  function  for  a wlde-band 
silicon  detector  (Fig.  9)  and  integrated  over  wavelength.  The  wavelength 
intervals  and  the  average  response  within  the  intervals  are  listed  in 
Table  III.  The  RADFLO  source  term  weighted  with  the  response  function 
and  integrated  over  wavelength  is  plotted  versus  time  out  to  24X10  ^ 
seconds  in  Fig.  10.  Figure  10  shows  the  response  that  should  be  observed 
if  atmospheric  attenuation  and  scattering  did  not  perturb  the  optical 
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Fig.  8.  Variation  of  the  Total  Flux  at  Satellite  with  Look 
Angle:  Ground  Albedo  =0.8 
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TABLE  III.  WIDE-BAND  SILICON  RESPONSE  FUNCTION 
VERSUS  WAVELENGTH 


WAVELENGTH  INTERVAL 
(micron) 

AVERAGE  RESPONSE 

0. 4000-0. A500 

0.1312 

0.4500-0.4850 

0 . 3000 

0.4850-0.5253 

0.4000 

0.5253-<^.5700 

0.4937 

0.5700-0.6266 

0.5875 

0.6266-0.6812 

0.6779 

0.6812-0.7358 

0.7548 

0.7358-0.8029 

0.8354 

0.8029-0.8700 

0.9125 

0.8700-0.9309 

0.9725 

0.9309-0.9918 

0.9800 

0.9918-1.0200 

0.8375 

1.0200-1.0600 

0.4820 

1.0600-1.1000 

0.1980 

[• 
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signal.  Note  that  first  maximum  occurs  at  approximately  7.08X10  ^sec- 
onds. The  optical  signal  seen  at  a receiver  position  Includes  not  only 
the  direct  but  also  the  scattered  radiation  and  since  the  scattered  radia- 
tion Is  delayed  In  time,  the  time  distribution  of  the  optical  signal 
at  the  receiver  position  may  be  quite  different  from  the  time  distribu- 
tion of  the  radiation  emitted  from  the  source.  Since  the  source  angular 
distribution  affects  the  ratio  of  the  scattered  to  direct  radiation 
observed  at  a receiver  position,  it  may  also  affect  the  time  distribution 
of  the  optical  signal  at  the  receiver  position.  That  is,  receivers  at 
the  same  look  angle  and  radial  distances  relative  to  the  source  but  at 
different  azimuthal  positions  will  show  a different  response  if  the 
source  is  anisotropic. 

Figures  11  and  12  show  the  calculated  time  response  for  receivers 
located  at  30-  and  90-degree  look  angles  from  each  of  three  sources  con- 
sidered. The  upper  curves  are  for  a source  emitting  two-thirds  of  Its 
radiation  in  the  forward  hemisphere  and  one-third  in  the  backward  hemi- 
sphere (FORWARD),  the  center  curves  are  for  an  Isotropic  source  (ISOTROPIC), 
and  the  lower  curves  are  for  a source  emitting  one-third  of  its  radiation 
in  the  forward  hemisphere  and  two- thirds  in  the  backward  hemisphere 
(BACKWARD) . Note  that  not  only  does  the  anisotropy  of  the  source  affect 
the  magnitude  of  the  optical  signal  received  at  the  detector  position, 
but  the  time  at  which  first  maximum  is  reached  also  appears  to  be  shifted 
to  later  times  as  the  source  becomes  peaked  in  the  backward  direction. 

The  shift  in  the  time  of  first  maximum  raises  a question  concerning  the 
time  scale  that  should  be  used  in  comparing  optical  signals  from  detectors 
located  at  different  receiver  positions. 

The  time  Intervals  used  in  the  POLY  program  to  record  the  time- 
dependent  optical  signal  are  in  the  order  of  10  to  20  microseconds  in 
w^dth  for  time  periods  near  the  time  of  first  maximum.  Therefore,  it 
is  difficult  to  determine  within  10  to  20  microseconds  when  first  maximum 
actually  occurs.  An  examination  of  the  optical  signals  generated  with 
the  POLY  program  indicates  that  scattering  may  shift  the  time  of  first 
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maximum  by  several  10 's  of  microseconds.  This  is  especially  true  if  the 
ground-reflected  signal  reaches  the  receiver  position  after  the  peak  in 
the  direct  radiation. 

In  order  to  determine  the  time  scale  that  should  be  used  in 
comparing  optical  signals  from  different  receivers,  the  time  of  1/2 
maximum  power  and  1/10  maximum  power  were  determined  from  the  POLY  data. 
These  times  are  shown  in  Tables  IV  and  IX  for  ground  albedos  of  0.0  and 

0.1. 

The  time  at  which  the  average  power  occurs  in  the  ith  time 
interval  varies  with  the  time  interval.  Since  it  was  not  possible  to 
obtain  the  time  at  which  the  average  power  occurred  in  each  of  the 
time  intervals  directly  from  the  POLY  calculations,  it  was  felt  that 
it  was  better  to  assume  for, the  purposes  of  this  study  that  the 
average  power  in  the  1th  time  interval  occurred  at  a time 


*^1-1 


(t^-t^_l)S,n2 


for  times  before  first  maximum  instead  of  at  a time  given  by  the 
midpoint  of  the  time  interval.  The  time  of  1/2  maximum  power  is  then 
given  by  logarltlimlc  interpolation: 


where  and  t^  are  the  times  of  the  two  powers  P(tj) 

bounding  the  power  P]^/2max  ~ maximum  power.  The  time  of  1/10  maximum 
power,  is  computed  by  replacing  the  power  at  1/2  max,  P2^^2max’ 

in  the  above  equation  with  the  power  at  1/10  maximum  power 
The  times  of  1/2  maximum  power  and  1/10  maximum  power  have  been  determined 
for  the  unattenuated  source  power  time  distribution  shown  in  Fig.  10. 

The  time  of  1/10  maximum  power  was  determined  to  be  39.09  microseconds 


30 


TABLE  IV.  POWER  AT  FIRST  MAXIMUM,  TIMES  OF  ONE-HALF  AND  ONE-TENTH  FIRST 
MAXIMUM:  25  KM  METEOROLOGICAL  RANGE  ATMOSPHERE,  HEIGHT  OF 

SOURCE  - 0.5  KM,  GROUND  ALBEDOS  = 0.0  AND  0.1 


GROl'NC 

A L 3 F 0 C # 

0 . n 

Ground 

albedo  0 

0 . 1 

SCU«''F 

FIRST* 

T I MR  Of  ** 

TIM£  OF** 

FIRST* 

TIME  OF** 

TIME  OF** 

TYPE 

MAX 

1/2  yAX 

1 / 1 0 MAX 

MAX 

1/2  MAX 

1/10  MAX 

l.00*< 

amglC 

0 

0 degrees 

forward 

7S0 . 0 

5^  . 1 1 

40  .67 

666.0 

55.69 

41.09 

ISOTROPIC 

6Q6.fi 

66 . 28 

4C  . 70 

716. 0 

55 .96 

41.34 

pac^<waro 

469. 0 

56.68 

40.96 

F66 . 0 

56.38 

41.73 

lock 

A Aj  U L E 

U 

30  DEGREES 

F0RU.AHD 

72P  . ? 

65.27 

4 0.73 

P14  . 3 

55.72 

41.15 

ISOTROPIC 

663 . 9 

65.43 

40.86 

647,6 

55.99 

41.38 

PAC<WARn 

399 . 6 

55 . 72 

41,09 

480.6 

5 6.44 

41.77 

look 

ANGLF 

» 

60  DEG9EES 

FORWARD 

639 . 3 

56 . 69 

41,23 

67  7.7 

55.91 

41.45 

I SOTRCP  I r- 

412.fi 

55.81 

41.3? 

449.1 

56.10 

41.60 

BACKWARD 

266 . 2 

56.03 

41  ,50 

320.6 

66  . 45 

41.89 

look 

A M G L E 

it 

70  DEGREES 

forward 

42  7.2 

55 .95 

41.58 

4 48.1 

56.11 

41.72 

I SOTROP  I D 

326 . 6 

56.07 

41.66 

346.0 

56.26 

41.85 

SACKWARH 

226.9 

56 . 27 

41.83 

243.9 

56 . 54 

42.08 

lock 

amgle 

it 

60  ■'ESREES 

FORWARD 

238.7 

56.47 

4?  .32 

245.6 

56.56 

42.41 

I SOTRDP I c 

181  . 9 

56.55 

42.42 

133.2 

66.72 

42.51 

R ACKWARD 

126  . 1 

56.82 

4?  .54 

1 3l  .2 

57. 26 

42.7? 

LOOK 

AMGL  F 

It 

90  DEGREES 

forward 

in  .31 

72.35 

62.31 

10.84 

72.  87 

52.59 

ISOTROPIC 

5 .229 

73.39 

52.77 

8 . 725 

74  . r 1 

53.11 

BACKWARD 

ft.  192 

76.27 

53.67 

6.660 

77.56 

54.10 

2 

* (watts/cm  ) 

**  (usee) 
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TABLE  V.  POWER  AT  FIRST  MAXIMUM,  TIMES  OF  ONE-HALF  AND  ONE-TENTH  FIRST 
MAXIMUM:  25  KM  METEOROLOGICAL  RANGE  ATMOSPHERE,  HEIGHT  OF 

SOURCE  - 3.0  KM,  GROUND  ALBEDOS  - 0.0  AND  0.1 


ground 

ALHEDP  tt 

0.0 

GROUND 

albedo  a 

0 . 1 

SOURCE 

FIRST* 

TIME  OF** 

T I a'  E OF  ** 

F IRST* 

TIME  OF** 

time  or** 

type 

AX 

1/2  MAX 

1/1 0 MAX 

MAX 

1/2  MAX 

1/10  MAX 

LOOK 

AriGLF 

» 

0 DEGREES 

rORRARH 

79b  . 5 

5a  . 94 

40.35 

855 . 0 

55.66 

40.85 

I SOTRORT  C 

6^3  . P 

55. 1? 

40  , 4R 

659.7 

56.17 

41.20 

BACKI-.  ARH 

Abl  . 1 

55.46 

40.73 

529.6 

57.67 

41.88 

LOOK 

A M G L E 

it 

30  DEGREES 

FORUARi) 

74C  . 4 

55.15 

4 0.56 

7 9 0.5 

55.81 

41.01 

ISOTROPIC 

R 7 n . I 

55 . 32 

40.68 

622 . 2 

56.20 

41.30 

B ACKWAwr; 

400.6 

55 .65 

40.94 

46^'  .1 

57.57 

41.93 

LOOK 

A i . G L E 

tt 

60  DEGRPES 

forrarp 

590.6 

55,67 

41,04 

612.8 

55.99 

41.27 

I soTROP [ r 

45l  . a 

55.50 

41.14 

474.5 

56.25 

41.46 

R ACKU  AR!) 

'5l4. 3 

56.10 

41  .37 

3^5.3 

56.97 

41.56 

LOOK 

angle 

p 

70  DEGREES 

FOR  ■JARP 

503.9 

55.91 

41  .34 

516.7 

56. 11 

41.49 

ISOTROPIC 

3B4 . 5 

56.03 

41  .43 

397.1 

56. 29 

41.63 

BACKRARP 

P65.9 

56 .28 

41.63 

279.4 

56.82 

41 . 95 

LOOK 

A'.GLE 

n 

80  degrees 

FOR'aARIj 

34  4.3 

56 .28 

4?  . 0 0 

349.0 

56 . 37 

42.07 

ISOTROPIC 

9 61.9 

56.37 

42.07 

266 . 2 

56.49 

42.16 

RACKW ARR 

179.7 

56.56 

42.22 

15  4.2 

56.90 

42 . 37 

LOOK 

angle 

« 

90  DEGREES 

F ORRARD 

39  . a 1 

60.03 

44 .62 

4 n C . 5 

60 .12 

44.66 

ISOTROPIC 

n . 0? 

60 . m 

44.66 

302.5 

60 . 21 

44.71 

R ACKWARC 

^0  . 3? 

60.31 

44.76 

206 . 2 

60.54 

44.86 

2 

* (watts/cm  ) 

**  (ysec) 
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TABLE  VI.  POWER  AT  FIRST  MAXIMUM,  TIMES  OF  ONE-HALF  AND  ONE-TENTH  FIRST 
MAXIMUM:  25  KM  METEOROLOGICAL  RANGE  ATMOSPHERE,  HEIGHT  OF 


SOURCE  = 

10.0  KM, 

GROUND  ALBEDOS 

= 0.0  AND  0.1 

GR  oor'n 

A|_3Er)0 

tf  0.0 

GRCi.iI  0 

ALEEdC  4 0 

, 1 

SOU'^C'? 

FIRST* 

time  of 

**  X I f)  r 

** 

r 1 RST  * 

TIK'E  CF** 

ll^ir-  Of-** 

TYPE 

^ AX 

1/2  max 

l/io  :;ax 

••■AX 

1/2  MAX 

: /lo  max 

look 

angle 

» 

0 oegpfrs 

FOR'aARO 

902.4 

54.57 

40.03 

»02  .6 

54.57 

40  . '>3 

ISOTROPIC 

613.6 

54.65 

4P  . 09 

613.7 

54.65 

40.09 

BACh:w  ARP 

425  . P 

54.41 

40 .22 

425.2 

54. G2 

40.22 

L00»' 

angle 

A 

30  DEGREES 

FORinARO 

716.7 

54.64 

4 0.10 

747.0 

54.64 

40.10 

1 SOTROP  T C 

566.5 

54.71 

40 .15 

566.7 

54.71 

40.15 

RACKwARtj 

3 4 7.1 

54  . «6 

40.26 

387.4 

54.87 

40.27 

LOOK 

AmGLF 

# 

60  degrees 

rOR'-.ARD 

610.7 

54 .93 

40 .33 

613.0 

54 . 97 

40.36 

ISOTROPIC 

162.7 

55. 01 

4 0.39 

469.4 

55.17 

40 .50 

8ACKWAR  ,J 

316.2 

55.21 

40.53 

328 . 1 

55 .62 

40.81 

LOCK 

' angle 

70  DEGREES 

TOR  W A R 0 

554.5 

55.24 

40.59 

558.9 

55.33 

40.64 

I SOTROP  r r 

420.6 

55.33 

40.65 

424.4 

55.43 

40.72 

BACKvJAPO 

2‘47 . 2 

55 .5Z 

40.79 

293.2 

55.75 

40.95 

LOOK 

Angle 

# 

GO  DEGREES 

FOHlvARC 

149  . R 

55.51 

4H  . 9g 

451.2 

55.64 

4i.no 

ISOTROPI'- 

3 4 0 . 3 

55.67 

41.03 

311.4 

55.70 

41.05 

R A C K W A R 0 

231  . C 

55.80 

41.14 

2 32.1 

55.85 

41.17 

LOOK 

AL'GlE 

90  DEGRFRS 

F 0 R A R 

10  7.0 

56 .50 

42 . 59 

1 0 T . 0 

56.50 

42.89 

ISOTROPIC 

. 52 

56 .53 

4 2.61 

8 0 . 5 '-> 

56. 53 

4 2 . A 1 

RAC*\K  AR’O 

54  .P5 

56  .5« 

42.64 

54.07 

56.58 

4 2.65 

2 

* (watts/cm  ) 
**  (usee) 
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TABLE  VII.  POWER  AT  FIRST  MAXIMUM,  TIMES  OF  ONE-HALF  AND  ONE-TENTH  FIRST 
MAXIMUM:  3 KM  METEOROLOGICAL  RANGE  ATMOSPHERE,  HEIC3IT  OF 

SOURCE  = 0.5  KM,  GROUND  ALBEDOS  =0.0  AND  O.i 


GRr.tiNP 

ALHtrPO  tr 

0 . n 

GPnii>'0 

ALPEDO  » 

0.1 

?ni.ipr  F 

r I p ST  * 

Tp'C  CF** 

X j y p nr** 

FIRST* 

Tiyp  OF** 

TIME  OF*F 

T YK'" 

■-I  AX 

X/2  '^AX 

1/1 C 14X 

’AX 

1/2  PAX 

1/10  »^AX 

look 

alGl"^ 

0 oFGRt-es 

7S  7 . 5 

58 . 92 

41.39 

851.9 

56.41 

41.84 

ISOMJrjpjp 

A 7 3 . 1 

5 6 . 1 r 

4 1.55 

766 . 6 

56 .62 

42.04 

gACYVApn 

“>  1 8 . 6 

56 .32 

41.76 

6P1 .5 

57 . 16 

42.30 

I.OPK 

AMGl? 

n 

30  OEGPFES 

rOR..APD 

791.8 

55.96 

4 1 .49 

8 7 8.6 

56 .34 

41.86 

1 SCTRnp  I r 

A n . ] 

56 .21 

4 1.70 

7 1,1 . 5 

56.64 

42.12 

p A c -<  '-i  4 n 

A 7 4 . A 

56,63 

42.07 

546.0 

57 . 70 

42.58 

tr  OK 

A 0 G L F 

« 

Af.  OEuRFFS 

TCP'/  t-  ' 1 

73  8 . 8 

5') . 1? 

41.90 

5 6 6.8 

56.30 

42.11 

I S(jT.;np  T ' 

4??  . 4 

‘>6  . <5 

42.13 

446.9 

56.57 

42.37 

ts  AC'</  ‘ • 1 

K7.  ^ 

■ ‘ . 9 9 

4p  . 49 

329  . 

57.60 

42. 76 

Lorjp 

AMGlE 

u 

7 0 “ F 9 ,<  r p p 

F OR^Am 

' T >■  . 

’0.2? 

41.92 

392  . " 

56  . ,38 

12.0a 

I SO  I '• 

T 

\ 6.5’ 

42.14 

31.1.3 

56.75 

42.31 

PACK..  AP 

R 1 8 . 1 

' T . 5 2 

42.59 

?36 . n 

56.  02 

42.  Si 

1.  n 0 K 

A'  GL'"- 

a 

80  DFGFFFS 

r pp  A R' 

' \ 

60 .73 

44.14 

117.1 

60.09 

4 4.71 

I SOTPrp  : ;• 

')  ■ . UP 

*-0  . 6? 

44.53 

94 , 62 

61.07 

44,74 

e A c '-j  A H 

hy  .(  b 

>P  .35 

45.31 

7 2.87 

6 2.89 

4 S . 5n 

L nr  K 

A^  GLF 

n 

90  2EGPFR3 

FOh  AP' 

0 . 4 3 3 

79.81 

54.31 

5.942 

79.92 

58  . A.8 

iscvn.Ptc 

pp . 77 

59 . 7« 

4.259 

82.96 

6 0.12 

(J  Api'U.  AP'' 

4 . ? 1 n 

F 7.29 

62  .09 

4.605 

8 7.53 

62.39 

2 

* (watts/cm  ) 
**  (usee) 
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tabu:  VIII.  POWER  AT  FIRST  MAXIMUM,  TIMES  OF  ONE-HALF  AND  ONE-TENTH  FIRST 


MAXIMUM: 
SOURCE  = 

3 KM  METEOROLOGICAL  RANGE  ATMOSPHERE 
3.0  KM,  GROUND  ALBEDOS  = 0.0  AND  0.1 

, HEIGHT  OF 

CFOijf^r 

A L o c D 0 « 

0 . 0 

GROliFD 

AUBFCG  n 0 

. 1 

SCU»(T: 
T vpc 

F ' HS  1* 
•AAX 

T 1 Qp** 

\ / ? "AX 

TIME  OF** 
l/lO  MAX 

f I R S T* 
'■AX 

TIME  CF** 
1/2  MAX 

riMc  OF** 

1/10  max 

LOOK 

AMGLF 

U 

0 5 R G 6 h 1-  S 

1 SOTPCP ! r 
S A c K w A P 

■<0.6 
‘ T i . 1 
o;>0  . -R 

53.41 
55.9’ 
66  . iiH 

40  .73 
41.05 
4 1.66 

6-i7 . 4 
7 3 ’ . 2 
5^2.3 

5 6.11 
56 . 90 
59.35 

41.17 

41.63 

42.47 

LO  O'' 

A >tGL F 

H 

30  nPG-'’RRS 

F ORWA^n 

ISOTROPIC 

gACKwARO 

796  . g 
6 A 0 . r 
493.9 

6=.  .07 
66.63 

59.05 

4 1.46 
41.56 
4 2.39 

641.0 
6 9 T . 0 
545.4 

56,51 
58 . 0« 
61.10 

41.52 

42.09 

43.18 

look 

AmGlE 

It 

60  DEGREES 

F OKWA  R^ 
ISOTROPIC 

R A c K Fi  A 0 [) 

64f>  .b 
^■>6 . 7 
370.4 

5b  .41 
56.91 
58.56 

41.68 

41.93 

42.50 

6 64.5 
52482 
389.9 

56. 76 
57.67 
59.55 

41.84 
42.16 

42.85 

LOOK 

AMGl  F 

it 

70  DEGREES 

F DP  V A P 
I SOTFJOP  1 : 
HACK  WAR'.' 

5o2  , 

136.7 

328.1 

56.50 
57.31 
59 . 48 

41  . fl4 
42.17 
42.95 

560.4 
4 45.4 

337.5 

56,64 
57 . 69 
59.99 

41.93 
42.30 
43. 14 

LOOK 

amgle 

tf 

80  DEGREES 

F OR  WARP 
i SOTPcR 1 C 
R A CK  A '' 

368.1 

281.6 

196.6 

56.42 

56.56 

57.1? 

42.17 
42.27 
42 .53 

370.1 
284 . 0 
199.3 

56.47 
56 . 63 
57.43 

42.20 
42 . 33 
42.6? 

LOOK 

A-’GLF- 

tt 

90  DEGREES 

F ORaiAR’; 

isotropic 

3ACKWAP0 

16  . 
11.96 
' . 9=1 

66.54 

67.63 

70.18 

48 . 66 
49.12 
50.06 

15 . 20 
12.16 
0.174 

66.69 

68.03 

70.61 

43.74 

49.25 

50 . 25 

2 

* (watts/cm  ) 

**  (yaec) 
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TABLE  IX.  POWER  AT  FIRST  MAXIMUM,  TIMES  OF  ONE-HALF  AND  ONE-TENTH  FIRST 
MAXIMUM:  3 KM  METEOROLOGICAL  RANGE  ATMOSPHERE,  HEIGHT  OF 

SOURCE  - 10.0  KM,  GROUND  ALBEDOS  - 0.0  AND  0.1 


GROUND 

ALSrDO  i, 

P . 0 

GRUUMH 

ALDEDO  a 

0 , 1 

SOUKct 

first  * 

TiMp  OF** 

Tjme  OF 

Irk 

FIRST* 

TIME  OF** 

TIME  Cf*^ 

type 

MAX 

1/2  MAX 

1/10  max 

MAX 

1/2  MAX 

1/10  MAX 

LOOK 

angle 

# 

0 DEGREES 

FCHMARn 

794.5 

54 .56 

40.03 

794 .6 

54,56 

40.  18 

I SQTROP IC 

606.1 

54.63 

40.08 

619.2 

54.87 

40 . 24 

HACKW APD 

429,  fi 

55.05 

4C  . 4 0 

474.6 

56.19 

41.16 

LOOK 

angle 

» 

30  DEGREES 

F 0RW4  Rn 

733.9 

54.38  . 

39,91 

785 . 4 

55.15 

40.43 

ISOTROPIC 

57B.6 

54,90 

40.26 

633.1 

55.91 

40 . 95 

BACKwART 

436.5 

56.16 

41.15 

486.1 

58.42 

41.99 

LOOK 

A N G L E 

# 

60  DEGREES 

F ORUARP 

633.3 

55.26 

40.57 

6 4 3.7 

55. 49 

40.70 

ISOTROPIC 

454.9 

55.50 

40  . 7l 

495.1 

55.73 

40.87 

P ACKWAPO 

33  6.9 

55 .96 

41.03 

349.3 

56.29 

41.26 

LOOK 

angle 

u 

70  DEGREES 

F ORWA  RD 

573  . ? 

56.10 

41.13 

598.5 

56.20 

41.25 

I SOTPOP  I C 

•15?  . 

56.27 

41  .30 

4 S 3 . 1 

56.39 

41.39 

BACKhArO 

314  . 1 

66.65 

n . 58 

319.6 

57. 08 

41.71 

LOOK 

A ' GL  E 

« 

8P  DEGPEFS 

T OR/JARr' 

472,8 

55.99 

11  .23 

474.1 

5 6,02 

41.25 

I SOTPnp  1 <- 

3S7.  B 

56.06 

41  .28 

5 5'-  .6 

56  . OV 

4l  . 3P 

R A C K W A R ■ 

:’<!?. 9 

56.70 

4l  .38 

2 M . 9 

56.  ?4 

4 1.41 

LOOK 

A^'^I|  r 

u 

OP  DEGorre; 

F 0 R P A P 

10  8,4 

56 . 35 

42 .3« 

n f- . A 

56.35 

42,38 

I SOTPOP  I c 

.53 

56.36 

4 2.30 

81.54 

56.37 

42.39 

b AnYwAPC) 

54 ,6V 

56.41 

12.43 

54,69 

56.41 

42. 43 

* (watts/cm^) 

**  (usee) 
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and  the  time  for  i/2  maximum  power  was  determined  to  be  51.72  microseconds. 
It  can  be  seen  in  Tables  IV  and  IX  that  for  those  cases  where  the  direct 
component  of  the  radiation  is  dominant,  l.e.,  small  look,  angles  and  high 
source  altitudes,  the  time  of  1/2  maximum  power  is  approximately  55 
microseconds  and  the  time  of  1/10  maximum  power  is  approximately  40 
microseconds.  The  times  to  1/2  maximum  power  and  to  1/10  maximum  power 
tend  to  increase  in  those  cases  where  scattering  predominates. 

The  source  angular  distribution  has  a small  effect  on  the  times 
of  1/2  maximum  and  1/10  maximum  power.  For  small  look  angles,  the 
peaking  of  the  source  in  the  backward  direction  shifts  the  time  of 
1/2  first  maximum  by  approximately  0.85  microseconds  and  the  time  of 
1/10  first  maximum  by  approximately  0.52  microseconds.  For  the  90- 
degree  look  angle,  the  average  shift  in  the  time  of  1/2  first  maximum 
is  2.66  microseconds  and  the  average  shift  in  the  time  of  1/10  first 
maximum  is  1.13  microseconds;  however,  for  those  cases  where  the  source 
is  at  500  meters  altitude,  the  shift  in  time  of  1/2  first  maximum  is 
as  much  as  7.48  microseconds  and  the  shift  in  the  time  of  1/10  first 
maximum  is  as  much  as  3.78  microseconds. 
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V.  POLO  CALCULATIONS  FOR  A 40-KM  METEOROLOGICAL 
RANGE  ATMOSPHERE 

POLO  calculations  have  been  previously  reported  (Ref.  14)  for 

receivers  at  synchronous  satellite  positions  and  sources  within  model 

atmospheres  having  ground  meteorological  ranges  of  3,  10,  and  25  km. 

During  the  reporting  period  additional  POLO  computations  have  been 

made  for  a model  atmosphere  with  a 40-km  meteorological  range  at  ground 

level.  The  POLO  calculations  were  performed  for  source  altitudes  of 

0.1,  0.5,  3,  and  10  km,  for  source  wavelengths  of  0.4278,  0.5,  0.6, 

0.75,  and  1.07  pm,  and  for  receiver  look  angles  of  0,  30,  45,  60,  70, 

75,  80,  and  90  degrees.  Ten  ground  albedo  values  ranging  from  0.0  to 

0.9  were  used.  The  POLO-calculated  scattered  fluxes  were  computed  for 

—8  —2 

delayed  times  between  10  and  10  secs.  The  source  was  taken  in  POLO 
to  be  a point-isotropic  instantaneous  emitter.  A description  of  the 
atmospheric  model  and  of  the  information  made  available  in  this  new 
data  set,  as  well  as  a discussion  of  the  usefulness  of  the  data,  are 
contained  in  this  section. 

5.1  Choice  of  Wavelengths 

Wavelengths  had  to  be  chosen  such  that  an  integration  over 
wavelength  for  a realistic  source  distribution  and  receiver  response 
in  the  TMTAU  procedure  would  give  accurate  data  without  use  of  excessive 
computer  time.  Using  data  computed  by  use  of  the  RADFLO  computer  pro- 
gram for  a low-yield  device  and  assuming  a detector  response  as  given 
in  Fig.  9,  the  product  of  the  RADFLO  data  and  the  receiver  response 
function  was  calculated  as  a function  of  wavelength  for  different 
emission  times.  The  resulting  data  are  plotted  in  Fig.  13.  It  can  be 
seen  from  Fig.  13  that  the  relative  effectiveness  of  the  radiation 
reaching  the  detector  system  varies  greatly  with  both  wavelength  and 
time.  For  early  time  Intervals,  most  of  the  radiation  being  detected 
is  in  the  range  between  0.4  and  0.8  pm  with  the  maximum  in  the  signal 
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occurring  for  wavelengths  in  the  0.44  ym  to  0.49  ym  interval.  For 
later  times,  the  wavelength  at  which  the  maximum  signal  occurs  will 


shift  to  longer  wavelengths.  For  larger  yield  devices,  it  is  expected 
that  the  wavelength  at  which  the  maximum  signal  occurs  for  a given 
time  will  shift  to  longer  wavelengths  since  the  smog  absorption  at 
small  wavelengths  is  expected  to  increase  with  increasing  yield. 

It  is  expected  that  the  POLO-calculated  transmission  data  would 
vary  smoothly  with  wavelength.  The  decrease  in  the  scattering  and 
extinction  cross  sections  result  in  scattered  transmission  data  that 
reduce  in  magnitude  with  increasing  wavelength,  whereas  the  direct 
transmission  tends  to  increase  in  magnitude  with  increasing  wavelength. 
Tae  magnitude  of  the  total  transmission  cata  changes  more  rapidly  at 
shorter  wavelengths  than  it  does  for  the  longci  wavelengths.  Therefore, 
it  was  decided  that  the  POLO  problems  for  the  40-km  meteorological 
range  atmosphere  should  be  run  for  wavelengths  of  0.4278,  0.5,  0.6, 

0.75,  and  1.07  ym. 

5.2  Model  Atmosphere 

The  U.  S.  Standard  Atmosphere,  1962  (Ref.  15)  was  used  to  cal- 
culate the  molecular  scattering  coefficient  as  a function  of  altitude. 
The  Elterman  1968  (Ref.  16)  model  was  chosen  to  describe  the  aerosol 
content  of  the  atmosphere  above  5-km  altitude.  The  value  of  the  aerosol 
coefficient  at  sea  level  for  0.55  ym  was  chosen  to  give  a meteorological 
range  of  40  km.  The  dependency  of  the  aerosol  scattering  coefficient 
between  0.0  and  5-km  altitude  was  assumed  to  vary  exponentially.  The 
dependency  of  the  ozone  absorption  coefficient  upon  wavelength  and  alti- 
tude was  taken  from  Ref.  9. 

Deirmend j ian ' s Haze  C model  (Ref.  17)  was  assumed  to  describe 
the  aerosol  size  distribution  for  all  altitudes.  MIE  calculations  for 
this  size  distribution  have  been  performed  previously  (Ref.  13);  the 
dependency  of  the  aerosol  scattering  and  the  aerosol  al sorption  coeffi- 
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dents  for  0.5,  0.6,  and  0.75  pm  were,  therefore,  taken  from  Ref.  13. 

The  data  for  0.4278  [im  and  1.07  ym  were  calculated  from  Ref.  13  under 
the  assumption  that  the  coefficients  vary  exponentially  with  wavelengths 
between  0.4  ym  and  1.55  ym. 

The  shape  of  the  curve  giving  the  wavelength  dependence  of  the 
aerosol  scattering  coefficient  at  sea  level  for  the  40-km  meteorologi- 
cal range  atmosphere  was  found  to  be  slightly  different  from  that  used 
previously  (Ref.  14)  for  the  3-,  10-,  and  25-km  meteorological  range 
atmospheres.  This  difference  results  from  the  fact  that  for  the  40-km 
meteorological  range  atmosphere  the  variation  of  the  aerosol  scattering 
coefficient  with  wavelength  was  based  on  the  MIE-2  calculations  for  the 
Haze  C aerosol  size  distribution,  whereas  for  the  3-,  10-,  and  25-km 
meteorological  range  atmospheres  the  wavelength  dependence  of  the  aero- 
sol cross  section  was  taken  from  the  data  reported  by  Elterman  for  a 
"clear  standard  atmosphere"  (Ref.  16). 

The  angular  scattering  data  for  the  40-km  model  and  for  0.5, 

0.6,  and  0.7  ym  were  taken  from  new  MIE  calculations  (Ref.  13).  For 
0.4278  and  1.07  ym,  the  phase  functions  that  were  used  for  the  problems 
in  Ref.  14  were  used  for  the  calculations  for  the  40-km  model  also. 

Water  vapor  absorption  for  1.07  ym  was  assumed  to  be  negligible 
in  the  calculations  for  the  40-km  model. 

5.3  POLO  Problem  Parameters 

There  were  180  problems  run  with  the  TPOLOW  version  of  POLO. 

Four  source  heights  (0.1,  0.5,  3.0,  and  10.0  kilometers),  five  source 
wavelengths  (.4278,  .5,  .6,  .75,  and  1.07  microns)  and  eight  look  angles 
(0,  30,  45,  60,  70,  75,  80,  and  90  degrees)  constituted  a set  of  160 
problems  which  were  run  for  two  time  ranges.  One  hundred  sixty  POLO 
problems  were  run  for  the  first  time  range  while  only  20  problems  were 
run  for  the  second  time  range.  Each  POLO  problem  computed  the  direct 
intensity  and  time-dependent  scattered  intensities  for  ten  values  of 
the  ground  albedos  (0.0  to  0.9  in  increments  of  0.1). 


The  two  time  ranges,  each  containing  50  intervals,  were  chosen 

such  that  an  overlap  of  30  time  Intervals  occurred.  The  earlv-time 
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range  was  from  0 to  10  second  with  an  upper  bound  for  the  first  time 
interval  of  one  shake.  All  problems  run  in  this  time  range  were  run 
for  one  look  angle,  5,000  histories,  a maximum  of  20  collisions,  and 
five  ground  reflections. 

A minimum  weight  cutoff  of  5X10  ^ was  used  to  terminate  histories 
for  those  photons  whose  weight  was  reduced  below  the  cutoff  value  where 
the  photon  weight  is  defined  as  the  fraction  of  the  source  photon  that 
remains  after  undergoing  absorption  and  scattering  interactions.  Biasing 
was  used  with  the  source  angular  distribution  to  force  more  of  the  source 
photons  to  leave  the  point  source  along  directions  toward  the  receiver 
position  and  to  force  collisions  to  occur  within  an  ellipsoid  with  foci 
at  the  source  and  receiver  positions  so  that  the  photons  would  arrive 
at  the  receiver  within  a retarded  time  of  5X10  ^ sec.  Biasing  of  source 
directions  and  of  collision  locations  was  used  in  POLO  to  provide  the 
best  statistical  results  for  the  early-time  portion  of  the  time-dependent 
scattering  calculations. 

-2 

The  second  time  range  was  from  0 to  10  second  with  the  upper 
bound  of  the  first  interval  at  10  ^ second.  All  problems  run  in  this 
time  range  were  run  for  all  eight  look  angles,  1,500  histories,  a maximum 
of  20  collisions  and  five  ground  reflections,  a minimum  weight  of  5X10 
and  no  biasing.  Comparison  of  data  from  the  two  time  ranges  demonstrated 
good  agreement  throughout  the  overlap  region.  The  average  ratio  of  the 
early-time  group  to  the  late-time  group  in  the  last  overlap  time  inter- 
val which  contained  information  (l.e.,  4.5X10  ^ < t ^ 5X10  was  .996 
for  all  problems  with  a 0-degree  look  angle  and  .821  for  all  problems 
with  a 90-degree  look  angle.  Results  for  these  two  look  angles  are  the 
best  and  worst  (respectively)  statistically.  A combined  set  of  70  time 
intervals  was  developed  by  using  data  from  the  early-time  group  for  the 

first  43  time  intervals  (0  to  5X10  ^ second)  and  the  last  27  time  inter- 
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vals  of  the  late-time  group  (5X10  to  1X10  second). 
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An  eight-digit  problem  number  was  associated  with  each  of  the 
180  problems.  Problem  numbers  were  of  the  form  MMLLLNJJ  where  MM  was 
the  meteorological  range  in  kilometers  (always  40) , LLL  was  the  wave- 
length in  hundredths  of  a micron,  N indicated  the  source  height  (N  = 0, 
1,  2,  3 corresponding  to  0.1,  0.5,  3.0,  and  10,0  kilometer  altitudes) 
and  JJ  was  the  look  angle  in  degrees  (except  for  JJ  = 99  which  indicated 
all  eight  look  angles  were  contained  in  one  problem).  Problem  numbers 
of  data  formed  from  combining  on  time  were  the  same  as  the  early-time 
problems . 

The  direct  flux  and  scattered  flux  integrated  over  time  for 
albedos  of  0.0,  0.1,  0.2,  0.4,  and  0.8  are  tabulated  in  Tables  X,  XI, 
XII,  and  XIII  as  a function  of  look  angle  and  source  wavelength  for 
source  heights  of  0.1,  0.5,  3.0,  and  10  km,  respectively. 

5.4  Variations  with  Meteorological  Range 

Comparison  of  the  POLO  calculations  for  a 40-kilometer  meteoro- 
logical range  with  the  data  previously  generated  for  the  3-,  10- , and 
25-kllometer  meteorological  ranges  shows  for  source  wavelengths  of 
0.4278  and  0.75  dm  that  smooth  semilogarithmlc  curves  can  easily  be 
fitted  to  the  data  as  a function  of  meteorological  range.  The  data 
for  the  four  different  meteorological  ranges  can  be  used  to  generate 
other  meteorological  ranges  without  additional  POLO  calculations.  Plots 
of  this  data  are  given  in  Figs.  14  and  15,  Similar  plots  for  the  1.07- 
micron  data  are  given  in  Fig.  16. 

In  previous  POLO  calculations  for  meteorological  ranges  of  3, 

10,  and  25  km,  the  data  for  1.07  ym  were  represented  by  an  average  over 
the  wavelength  interval  from  1.025  ym  to  1.175  ym.  This  wavelength 
range  contains  the  H2O  (})-band  around  1.13  ym.  From  absorption  data 
given  in  Ref.  18,  it  was  estimated  for  a look  angle  of  0°  that  the 
4>-band  will  reduce  the  averaged  intensity  within  the  1.025  to  1.175  ym 
wavelength  Interval  by  approximately  a factor  of  2 over  that  which 
would  have  occurred  if  H2O  absorption  were  neglected,  which  Is  in  good 
agreement  with  the  data  shown  in  Fig.  16.  For  larger  look  angles,  only 
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Fig.  16.  Variation  of  POLO  Calculations  for  1.07-micron  Wavelength  Light 
with  Meteorological  Range 


small  additional  effects  of  the  H2O  band  are  expected  and  there  should 
be  no  appreciable  further  reduction  In  Intensity. 

The  data  shown  In  Ref.  18  reveal  further  that  there  Is  no 
sizeable  absorption  by  H^O  (or  any  other  absorbing  species)  In  the 
range  1.00  - 1.07  pm.  Even  though  the  previous  POLO  calculations  were 
accomplished  for  three  subintervals  ranging  from  1.025  pm  - 1.075  pm, 
1.075  p - 1.125  pm,  and  1.125  pm  - 1.175  pm,  no  printout  was  obtained 
for  the  data  within  these  small  subintervals  and,  therefore,  no  direct 
comparison  with  the  new  data,  which  neglected  water  vapor  absorption, 
can  be  obtained. 


5.5  Curve-Fit  Coefficients  for  Input  to  TMTAU 


To  use  the  POLO  data  for  a AO-km  meteorological  range  atmosphere 
In  the  TMTAU  code  requires  that  the  data  be  expressed  mathematically  as 
a series  of  exponential  terms.  Two  codes,  CFIT  and  DFIT,  were  developed 
to  curve  fit  the  POLO  data  In  this  manner.  The  fitted  data  Is  output 
In  a format  suitable  for  Input  to  TMTAU.  CFIT  and  DFIT  were  developed 
from  the  TPOFIT  code  which  was  used  previously  to  fit  the  POLO  data  for 
atmospheres  with  ground-level  meteorological  ranges  of  3,  10,  and  25  km. 


The  POLO  data  are  given  as,  Ig(t),  the  scattered  flux  at  the 
receiver  within  differential  time  bins.  The  curve-fit  codes  first  cal- 
culates the  cumulative  scattered  flux,  I(t),  as  a function  of  retarded 
time,  T,  by 


I(T) 


'T 

Is(t)dt 


0 


A typical  set  of  these  curves  Is  shown  In  Fig.  17.  The  family 

of  curves  shown  there  Is  for  a source  altitude  of  0.5  km,  a wavelength 

of  0.6  microns,  and  an  albedo  of  a » 0.9.  These  curves  were  plotted 

—8  —2 

from  70  data  points  for  times  ranging  from  10  to  10  seconds.  The 
receiver  angle  parameter  denotes  the  polar  position  angle  between  the 


PROBLEM  400602XXALBED0  0.9 
SOURCE  ALTITUDE  3.0  KM 
WAVELENGTH  0.60  MICRONS 
METEOROLOGICAL  RANGE  40  KM 


Example  Plot  of  the  Cumulative  Intensity  vs  Ketarded  Time  for  X - 0.6  micron 


radials  through  the  source  and  detector  positions.  A drastic  change 
in  slope  for  the  curves  plotted  in  Fig.  17  can  be  noted  at  the  point  of 
inflection,  for  the  curves  for  receiver  polar  position  angles  one 
through  seven.  It  is  at  this  time  that  the  reflected  Intensity  from 
the  ground  starts  arriving  at  the  detector.  The  dashed  line  in  Fig. 

17  shows  I(t)  for  receiver  polar  position  angle  number  one  (0^  = 0°)  and 
for  a ground  albedo  of  0.  In  fitting  the  POLO  data,  a maximum  of  six 
terms  were  used  for  the  zero  albedo  case.  A maximum  of  three  terms 
were  used  to  fit  the  reflected  contributions  past  time  T^.  Thus,  the 
data  are  expressed  as: 


I(T)  = K, 


i<6 


1=1 


-b  T 

aie  + 


j<3 

- I 

j=l 


a . e 
J 


(1) 


where: 

= maximum  I(t)  for  zero  albedo  case, 

K2  = maximum  I(t),  including  reflection,  minus 

= earliest  arrival  time  of  reflected  photons. 

It  is  the  function  of  the  curve-fit  codes  to  evaluate  the  a^,  a.,  b^, 
and  bj  constants  for  each  curve.  The  method  used  to  obtain  these  con- 
stants has  been  described  in  Ref.  3 as  a "partitioned  least  squares" 
method. 

Data  output  from  the  curve-fitting  codes  consists  of  two  files, 
print  and  punch.  The  print  file  is  usually  assigned  to  the  line  printer. 
A sample  print  file  output  is  shown  in  Table  XIV  and  listing  of  a sample 
set  of  the  punch  file  data  is  shown  in  Table  XV. 

The  sample  listing  in  Table  XIV  shows  only  the  first  35  of  the 
70  time  points  which  were  used.  Pairs  of  curve-fit  coefficients  are 
shovm  at  the  top  of  the  page.  Whenever  a pair  are  both  zero,  this  means 
that  less  than  the  maximum  allowable  number  of  terms  (6  for  zero  albedo. 
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TABLE  XIV.  SA.KPLE  LISTING  OF  THE  PRINT  FILE  FROM  THE  CFIT  CODE 
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9 for  other  albedo  cases)  were  used  to  fit  the  curves.  The  number  of 
terms  to  be  used  is  a user  input  option  to  the  curve-fitting  codes. 

The  column  headed  POLO  DATA  in  Table  XIV  is  the  cumulative  sum 
of  the  POLO  data;  the  column  FITTED  DATA  is  the  cumulative  sum  calcu- 
lated by  the  code  using  the  curve-fit  coefficients.  The  STANDARD  DEVIA- 
TION AND  PERCENT  ERROR  are  calculated  in  the  standard  manner  and  refer 
to  the  difference  in  the  POLO  DATA  and  the  FITTED  DATA  in  the  output 
from  the  CFIT  code.  If  the  output  is  from  the  DFIT  code,  then  they 
refer  to  the  difference  between  DIFF  DATA  and  DERIVATIVE.  The  DIFF  DATA 
is  calculated  numerically  from  the  POLO  DATA  as  AI/At,  whereas  DERIVATIVE 
data  is  the  analytic  derivative,  dl/dt,  of  the  curve  which  is  produced 
by  inserting  the  curve-fit  coefficients  in  the  above  equation.  Both 
codes  make  several  tries  at  generating  an  adequate  set  ol  curve-fit 
coefficients.  Each  code  selects  the  best  set  of  these  coefficients 
based  on  the  standard  deviation  calculated  over  all  time  points. 

The  problem  parameters  used  to  define  the  POLO  problems  for 
which  curve-fitted  data  were  generated  are  shown  in  Table  XVI.  The 
curve-fit  parameters  were  all  generated  using  the  CFIT  code.  An  eight- 
digit  problem  " mber  is  associated  with  each  of  the  160  problems.  Prob- 
lem numbers  are  of  the  form  MM  LLL  N XX  where: 

MM  is  the  meteorological  range  in  kilometers, 

LLL  is  the  wavelength  in  hundredths  of  microns, 

N is  associated  with  the  source  height  (N  = 0, 1,2,3 
corresponds  to  0.1,  0.5,  3.0,  and  10.0  km), 

XX  is  the  look  angle  in  degrees. 

The  POLO  problems  for  the  40-km  meteorological  range  atmosphere  were  all 
run  for  detectors  positioned  at  synchronous  satellite  altitude.  In 
Table  XVI  the  XX  number  always  varies  over  each  eight-problem  set  as 
0,  30,  45,  60,  70,  75,  80,  and  90.  For  each  problem  curve-fit  data 
were  obtained  for  ground  albedo  values  of  0.0,  0.2,  0.5,  and  0.8. 

The  punch  file  output  shown  in  Table  .XV  is  a sample  listing  of 
a selected  part  of  the  data  which  now  exists  on  punched  cards.  The 
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TABLE  XVI.  SEQUENCE  OF  CURVE-FITTED  POLO  DATA  FROM  CFIT  PUNCH  FILE 


40-KM  METEOROLOGICAL  RANGE 


SEQUENCE 

NUMBER 

PROBLEM 

NUMBER 

SOURCE 
ALT (km) 

WAVELENGTH 
OF  SOURCE (pm) 

1-8 

400430XX 

0,1 

0.43 

9-16 

400500XX 

0.1 

0.50 

17-24 

400600XX 

0.1 

0.60 

25-32 

400750XX 

0,1 

0.75 

33-40 

401070XX 

0.1 

1.07 

41-48 

400431XX 

0.5 

0.43 

49-56 

400501XX 

0.5 

0.50 

57-64 

400601XX 

0.5 

0.60 

65-72 

400751XX 

0.5 

0.75 

73-80 

401071XX 

0.5 

1.07 

81-88 

400432XX 

3.0 

0.43 

89-96 

400502XX 

3.0 

0.50 

97-104 

400602XX 

3.0 

0.60 

105-112 

400752XX 

3.0 

0.75 

113-120 

401072XX 

3.0 

1.07 

121-128 

4004 33XX 

10.0 

0.43 

129-136 

400503XX 

10.0 

0.50 

137-144 

40U603XX 

10.0 

0,60 

145-152 

400753XX 

10.0 

0.75 

153-160 

401073XX 

10.0 

1.07 

foraat  of  these  data  is  compatible  with  the  input  requirements  of  the 
TMTAU  program.  Note  that  in  the  problem  number,  the  number  AO  has 
been  dropped  as  a prefix.  The  angle  set  which  is  the  last  two  digits 
refers  to  the  look  angle  rather  than  receiver  angle. 

Curve-fit  coefficients  of  the  POLO  results  for  the  AO-kilometer 
visibility  atmosphere  were  generated  for  four  different  values  of  the 
ground  albedo,  0.0,  0.2,  0.5,  and  0.8,  and  these  data  now  exist  on 
tape  at  AFTAC.  Since  the  TMTAU  procedure  assumes  that  curve-fit  data 
for  ten  values  of  the  ground  albedo  are  on  tape,  the  remaining  data 
(data  for  ground  albedo  values  other  than  0.0,  0.2,  0.5,  and  0.8)  were 
given  values  of  zero  and  stored  on  the  tape  in  the  appropriate  places. 
Therefore,  the  existing  version  of  TMTAU  can  be  used  with  the  AO-km 
meteorological  range  data  if  the  ground  albedo  requested  is  one  of  the 
four  mentioned  above.  The  data  resides  on  tape  number  001706  (a  standard 
label,  TAPE9,  external  label  FORTYKM-RBL,  tape)  in  three  forms. 

Label  1 (DSN=F0RTY)  contains  card  Images  of  the  data  for  only 
four  values  of  the  ground  albedo.  Label  2 (DSN=FORTY10)  contains  card 
Images  of  the  data  spaced  for  all  ten  albedos.  That  is,  there  are  three 
blank  cards  Inserted  for  the  data  for  each  missing  value  of  the  ground 
albedo.  Both  label  1 and  label  2 have  DCB= (RECFM=FB ,LRECL=*80 ,BLKSIZE= 
3200) . 

Figure  18  shows  the  data  for  problem  A00A3300  as  it  resides 
under  label  1.  Label  3 (DSN-FORTYB)  contains  the  curve-fit  data  in 
binary  format  which  is  the  required  format  of  the  TMTAU  and  TMTRFOl  pro- 
cedures. Label  3 has  DCB-(RECFM-VS,LRECL=92 ,BLKSIZE»3200) . Problem 
numbers  for  the  data  are  those  discussed  above  with  the  exception  that 
the  card  image  data  shows  the  problem  number  modulo  AOOOOOOO . 

Figure  19  shows  a list  of  the  header  cards  from  the  160  prob- 
lems which  reside  on  tape  1706  with  the  addition  of  a sequencing  number. 
On  cards,  the  P in  PR03  is  in  column  one. 
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1 PflOB0433n0sREC 

2 PR  OB  04  33 3C RFC 

3 PROB043345REC 

4 PROB043360REC 
3 PROB043370REC 
o PRO0O43375REC 

7 PR03Q43380REC 

8 PR09043390REC 

9 PR03050300REC 

10  PR  OB  05  03  30  RE  C 

11  PROR050345REC 

12  PR0905036nREC 

13  PR  OR  05  03  70  RE  C 

14  PR  OB  05  03  75  RE  C 

15  PR  08  05  03  80  RE  C 

16  PROB050390REC 

17  PROB  06  03  00  REC 

18  PROB060330  REC 

19  PROB060345REC 

20  PROB060360REC 

21  PROB  06  03  70  REC 

22  PR09060375REC 

23  °ROB060360REC 

24  PR03a60390REC 

25  PR  OB  07  53  00  RE  C 

26  PR  OB  07  53  30  RE  C 

27  PR0B075345REC 

28  PR08075360REC 

29  PROB075370REC 

30  PROB075375REC 

31  PROB075380REC 

32  PROB075390REC 

33  PROB107300REC 

34  PROB107330REC 

35  PROB107345REC 

36  PROB107360REC 

37  PROB107370REC 
33  PR0B1C7375REC 

39  PR08107380REC 

40  PROBIO 7390REC 


1DIRI«5.5482E-11 
2D IR 1=5. 2199E-11 
3DIRI=4.8161E-ll 
4DIR 1=4. 2209E-11 
5DIRI=3.6232E-11 
6DIR 1=3. 1844E-11 
7D  IR  1 = 2 . 5356E-H 
80 IR 1=6. 5649E-13 
ID IR 1=5. 7141E-11 
2D  IR  1 = 5.  40  04E-11 
3DIRI=5. 0208E-11 
40  IR  1 = 4. 4766E-11 
5D1R  1 = 3.  94  76  E-31 
6D  IR  1=  3.  5650E-11 
7D  IP  1=  2.  9950E-11 
801R 1=2. 4613E-12 
ID  IR  1=  5. 6369E-11 
2D IR 1=5. 3709E-11 
3D IR 1=4. 9874E-i1 
4DIRI=4. 4352E-11 
5DIRI=  3.  8963E-11 
6D IR 1=3. 5069E-11 
7DIRI=2.9323E-31 
8D1R1=3.6454E-12 
ID  IR  1=  5.  99  63E-11 
2D IR 1=5. 7095E-11 
3D  IR  1 = 5. 37  47E-H 
4DIR1=4.9282E-11 
5DIRI=4.54G1E-11 
6DIRI=4.2332E-ll 
7DIRI=3.9178E-ll 
8DIRI=1.35l7E-ll 
1DIRI=6.0928E-11 
2D IR 1=5. 8157E-11 
3D IR 1=5. 4974E-11 
4DIRI=5. 0877E-11 
5DIRI=4. 7553E-11 
6PIRU4 .5516E-11 
7DIRI=4.2816E-11 
8D IR 1=2. 3025E-11 


IN 
M IN 
MIN 
M IN 
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M IN 
MIN 
M IN 
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M IN 
MIN 
MIN 
M IN 
M IN 
M IN 
M IN 
M IN 
M IN 
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IN 
M IN 

IN 
MIN 
MIN 
M IN 
M IN 
M IN 
M I N 
MIN 

IN 
MIN 
M I N 
M I N 
M IN 
MIN 
MIN 
MIN 
MIN 
MIN 
MIN 


D/S=  . 44flTCRIT  = l ,5000E*’i« 
n/8=  . 8n5TCR  iT=i  .5  00  0-:- -.-4 

D/8=  . 565'^CR  1 1 = 1 . 500CE- '18 

D/0=  .668TCRIT  = 1. 50003- iJ« 
D/S=  .829TC8lT=1.50COt-08 
D/S=  . ‘>57TCR  IT=1 .50003 8 
D/S=  . 98  2TCR 1 1=1 .5C0CE-O8 

D/8=  . P04TCR 11=1  .50C0E-0C 

D/S=  . 53  0'"CRI1=1  .5C'00E-0e 
D/S=  . 57  6T  OR  IT  =1  .5  CO  OE -0  8 
0/S=  . 652TCRH=1  .500  03-08 

P/S=  . P34TCR  11  =1  . 5 00  CE-Cfe 
D/?=  1.  04  2T  CR  II  =1  .5  CO  OF -0  e 
n/S=  1.  leOTCR  11  =1  .5  00  0E-08 
D/S=  1.  35  51  CRH  =1  .5  00  0b-0d 
D/S=  . P34TCR  11  =1  .5  00  0F -08 
D/S=  . 558TCR  11  =1  .5  00  0E-08 
D/'^=  . 625TCR  II  =1  .5  00  OE -06 

D/S:  . 7301  CR  1 1 =1  .5  00  Of -08 

C/8:  . 9431CR 11=1 .5COOc-On 

D/S=  1.  17  8TCR  II  =1  .'-■•.OOE-Oe 
D/S=  1 . 406TCR H =1 .50COF-"f 
D/S:  3 . 6901CR 1 1 =1 . 50C0F -06 
n/S=  1.  66  3TCR  11  =1 .5  00  00-0  8 
D/S:  . 56  4T  CR  II  =1  .5  CO  OF -C  b 

D/S:  . 64  ITCR  II  =1  .5  00  OE -0  8 

D/S:  . 757TCR  II  =1  .50C0E-0e 

D/S=  3 . 029TCO  1 1 =1  . 5000E-0  8 
D/S=  1. 3?6TCRI1=1  .5C00F-0o 
D/S=  3 . 600TCR 11  = 1 . 5C0CF -06 
D/S=  1. 9821CR 11=1 . 5C00E-0n 
D/S=  3.30STCRH  = 1.5CC0F-0i- 
D/S=  . 591 tcR 1 1=1. 50  OOF -Co 
D/S:  . 676TCR  11=1 .5''OCE-06 

D/S=  . 808TCR  II =1  .5  0CCE-06 
D/S=  1 . 106TCR  11  = 1 .5C0CE-0P 
D/S:  1 . 494TCR  1 1 = 1 . 5C00F -08 
D/S=  3.746TCRll=1.5CC0r-Co 
D/S=  2. ?64TCRI1=1 .50C0E-CP 
D/S=  4 . 699TCR  II =1  . 5CF OE-0  8 


Fig.  19.  Header  Cards  for  Curve-Fit  Data  on  Tape  1706 
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41  pRoBfl43100  REC  lDIRla3.92l6E-ll  'IIN  D/S=  . 2?  6T  CR  I T =1  . 5 CO  OE -0  6 

42  PROB043l3nREC  2DIR  1 =3  4978E-H  D/S=  2 4 3T  CR  I T = 1 . &C  qO  E- o8 

43  PROR043145REC  3D  I R I = 2 95  1 lE- 1 1 M j N D/''-  2^  9 T CR  i T = 1 . 5 0 0 OE  - 0 fi 

44  PROB043160REC  4 D I R I = 2 . 1 14 9E  - 1 1 MIN  D/S=  . 2P 7 T CR  I 7 = 1 . 5 0 0 0 E - 0 8 

45  PR05043170REC  5D I R I = 1 . 3254 E - 1 1 MIN  D/S=  . 2 7 0 T CR I T = 1 . 5 0 C 0 E - 0 8 

46  PROB043175REC  6D I R I = 3 . 4988E- 12  min  d/S=  . 24 4 T CR I T = 1 . 5 C C OE - 0 b 

47  PROB043180REC  7D IR I *3 . 6527E-12  MIN  D/S-  . 1 7 1 T CR  I T = 1 . 5 0 0 OE -0 8 

48  PROB04319OREC  8D I R I = 3 . 19 24 E- 20  MIN  D/S=  . 00 OT CR  I T =2  . 5 00  OE -0 6 

49  PROH050100REC  ID I R I = 4 . 49 2l£- 11  MIN  D/S=  . 28 1 T CR I T = 1 . 5 C 0 OE - L b 

50  PROB050130REC  2D  I R I =4 . 0915E-11  MIN  D/S=  . 31 4 T CR I T =1  . 5 0 0 OE -0 fc 

51  PROB050145REC  3D  I R I = 3 . 5754 E -11  MIN  D/S=  . 353TCR I 7 = 1 . 5000E-38 

52  PRDB050160REC  4D I R I =2 . 7730E-11  MIN  D/S=  . 4i i t CR 1 T = 1 . 50 0 OE - 0 « 

53  PROB050170REC  5D 1 R I = 1 . 9663E- l 1 MIN  D/S=  . 436T CR I T =1 . 5 0 0 OE -0 c 

54  PROB050175REC  6D I R I * 1 . 42 7 1 E - 1 1 MIN  D/S=  . 44  0 T CR I T = 1 . 5 0 0 OE  - 0 8 

55  PR08050130REC  7D IR I = 7 . 81 77 E- l2  MIN  D/S=  . 39 4T CR I T =1  , 5 00  OE -0 P 
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VI.  METHODS  FOR  FOLDING  TIME-  AND  SPECTRAL-DEPENDENT 
SOURCE  AND  ATMOSPHERIC  TRANSMISSION  DATA 


The  two  programs,  TMTAU  (Ref.  3)  and  POLY  (Ref.  1),  were  pre- 
viously developed  for  use  in  folding  a time-  and  wavelength-dependent 
source  term  with  atmospheric  transmission  data  generated  with  the  POLO 
program.  These  programs  have  both  been  improved  during  the  past  year. 
Both  programs  were  originally  developed  to  determine  the  time-dependent 
response  of  a wide-band  receiver  to  the  radiation  transmitted  through 
the  atmosphere  after  being  emitted  from  a time-  and  wavelength-dependent 
source.  The  total  radiation  I(x),  as  a function  of  time  T at  a receiver 
from  a source  S(X,t)  that  has  transmission  I(X,t),  and  a detector 
response  R(X) , where  X is  the  wavelength  and  t the  time,  is  calculated 
with  the  double  ..ntegral 


I(T) 


R(X)dX 


T 

S(X,t)I(X,T-t)dt 


X 


1 


0 


The  major  difference  in  the  T!ITAU  and  POLY  programs  is  that 
the  POLY  program  utilizes  the  source  and  transmission  data  expressed  in 
tabulated  form,  whereas  the  TMTAU  program  requires  that  the  transmission 
function  be  expressed  as  analytical  function  of  a particular  type. 

In  previous  versions  of  the  TMTAU  program,  the  transmission 
data  were  limited  to  that  for  three  wavelengths  and  the  time-dependent 
source  data  for  the  same  three  wavelengths  were  computed  by  TMTAU  from 
blackbody  scaling  laws  developed  by  T.  White  (Ref.  3).  Subsequently, 
the  TMTAU  program  was  modified  to  allow  the  use  of  other  source  data 
such  as  that  generated  with  the  RADFLO  program.  Later,  another  version 
of  TMTAU  was  developed  which  removed  the  limitation  that  the  source  and 
transmission  functions  be  defined  for  only  three  wavelengths.  This  new 
version  of  TMTAU,  designated  (TMTRFOl) , does  require,  however,  that  the 
source  and  transmission  functions  be  defined  for  the  exact  same  wave- 
length values. 
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Examination  of  the  wavelength  dependence  of  the  source  data 
generated  with  the  RADFLO  program,  which  is  currently  being  used  in  the 
TMTAU  and  POLY  programs,  lead  to  the  conclusion  that  the  source  data 
should  be  defined  at  more  than  three  wavelengths.  The  variation  with 
wavelength  of  the  atmospheric  transmission  data  generated  with  the  POLO 
program  is  much  smoother  than  the  source  variation  with  wavelength. 
Therefore,  a computer  program  designated  POLICE  was  written  to  generate 
the  curve-fit  coefficients  required  by  TMTRFOl  which  predict  time- 
dependent  transmission  data  for  wavelengths  other  than  those  obtained 
with  the  POLO  program.  A description  of  the  POLICE  program  is  given 
in  the  following  section  of  this  report. 

The  POLY  program  approximates  the  convolution  Integral  through 
applications  of  various  forms  of  numerical  integration  to  the  raw  POLO- 
calculated  atmospheric  transmission  data  and  a tabulation  of  the  product 
of  the  source  term  and  the  detector  response  function.  The  main  advan- 
tage of  the  POLY  program  is  that  it  eliminates  the  need  to  curve-fit 
the  atmospheric  transmission  data.  In  the  POLY  program,  the  times  and 
wavelengths  used  in  the  definition  of  the  source  term  do  not  have  to 
correspond  in  number  or  value  to  those  used  in  the  definition  of  the 
transmission  data.  Internal  to  the  POLY  program,  an  interpolation  of 
the  transmission  data  over  wavelength  is  performed  to  define  the  trans- 
mission versus  time  for  the  wavelength  values  used  in  the  definition  of 
the  source  terms.  A third  set  of  times  is  used  to  specify  the  particular 
time  values  for  which  the  time-dependent  power  at  the  detector  is  to  be 
determined . 

6.1  Code  for  Generating  Curve-Fit  Coefficients  for  Wavelengths 

not  Considered  in  POLO  Calculations 

The  POLICE  code  was  developed  to  generate  curve-fit  coefficients 
for  input  to  TMTRFOl  which  will  predict  time-dependent  scattering  data 
for  wavelengths  other  than  those  run  with  the  POLO  program.  The  POLICE 
code  aliminates  the  necessity  of  excessive  numbeis  of  POLO  runs  to  gain 
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a sufficiently  fine  wavelength  structure  for  the  data  used  in  TMTRFOl. 

The  code  is  variably  dimensioned  and  is  written  in  FORTRAN  for  use  on 
an  IBM  system. 

The  direct  flux  for  the  wavelengths  not  treated  by  POLO  is 
calculated  optionally  in  POLICE  by  interpolation  over  wavelength  of  the 
direct  fluxes  computed  with  POLO  or  by  use  of  the  subroutine  DBEAM. 

This  is  the  same  subroutine  that  is  currently  used  in  POLO  to  calculate 
the  direct  flux  at  a receiver  position. 

Curve-fit  coefficients  for  the  scattered  fluxes  at  wavelengths 
not  considered  in  the  original  POLO  calculations  are  optionally  calcu- 
lated in  two  ways.  Coefficients  used  to  define  the  "standard"  wave- 
lengths (0.4278  y,  0.75  y,  and  1.07  y)  are  interpolated  by  constructing 
an  interpolating  polynomial  of  degree  n-1  to  the  natural  logarithms  of 
the  coefficients,  where  n is  the  number  of  standard  wavelengths,  and 
then  evaluating  the  expression  at  the  desired  wavelengths.  The  inter- 
polating polynomial  is  constructed  using  the  divided  difference  method 
(Ref.  19).  Time-dependent  scattered  flux  data  for  the  standard  wave- 
lengths are  then  generated  from  the  coefficients  for  a time  base  which 
is  user-controlled.  Scattered  flux  data  for  the  standard  wavelengths 
are  then  interpolated  over  wavelength  for  each  time  in  the  same  fashion 
as  were  the  coefficients  and  the  results  temporarily  stored.  For  each 
new  wavelength,  a comparison  is  then  made  between  the  data  generated 
from  the  interpolated  coefficients  and  the  data  interpolated  from  the 
standard  data  evaluated  on  the  time  base.  If  the  comparison  is  acceptable, 
the  interpolated  coefficients  are  eventually  output;  if  unacceptable,  the 
data  which  were  interpolated  from  the  standard  data  evaluated  on  the  time 
base  will  be  retained  for  curve  fitting. 

An  option  in  POLICE  will  allow  the  user  to  calculate  the  scat- 
tered flux  at  wavelengths  other  than  X * 0.4278  y,  0.75  y,  and  1.07  y 
by  interpolating  on  a "build-up"  term  which  is  derived  from 

lg(t,X)  - lp(X)  B(t,X) 
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for  A • 0.4278  M,  0.75  u,  and  1.07  y,  where  t is  time,  A is  wavelength, 
Ig  is  the  scattered  flux,  1^  is  direct  flux  and  B is  the  build-up  term. 
The  scattered  fluxes  as  a function  of  time  are  generated  by  POLICE  for 
other  wavelengths  in  this  manner.  The  results  are  then  compared  with 
the  scattered  flux  data  generated  from  the  Interpolated  curve-fit 
coefficients.  If  the  comparison  is  acceptable,  then  the  interpolated 
curve-fit  coefficients  are  output.  If  the  comparison  is  not  acceptable, 
the  scattered  fluxes  generated  with  the  build-up  factors  are  retained 
for  curve  fitting. 

Three  separate  criteria  for  accepting  Interpolated  curve-fit 
coefficients  are  input  to  POLICE:  1)  the  minimum  ratio,  2)  the  maximum 

ratio,  and  3)  the  mean  ratio  of  the  data  when  considered  collectively 
over  the  time  base.  Should  any  criteria  not  be  met,  the  data  as  inter- 
polated from  the  standard  data  evaluated  on  the  time  base  will  be 
retained  for  curve  fitting  with  the  use  of  the  CFIT  or  DFIT  procedures. 


6.2  Checkout  of  the  TMTAU,  TMTRFOl  and  POLYN  Procedures 

An  exacting  checkout  of  the  TMTAU,  TMTRFOl  and  POLYN  procedures 
has  been  performed  to  insure  proper  operation  of  the  procedures  and  to 
establish  the  precision  of  the  results  obtained  from  the  three  computer 
procedures.  POLYN  is  an  improved  version  of  the  POLY  procedure. 

TMTRFOl  and  POLYN  have  been  established  at  AFTAC  on  the  disk  pack  TF2PAK 
under  RRACODES.  The  checkout  of  the  convolution  codes  was  performed  by 
comparing  program  results  with  analytic  solution  results.  The  compari- 
sons showed  good  agreement,  even  with  a nonreallstic  source  term. 

A fictitious  source  s(A,t)  was  defined  as 


2 

s(A,t)  = 3q  ^ e 


where  a^  = 3 and  a^  = 10  . A fictitious  transmission  term  I(A,x)  was 
defined  as 


I(A,x)  = QqA 


-a^^X 
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where  QIq  = 5 and  = 10^.  The  convolution  integral  is  defined  by 


F^(t)  F2(T-t)dt. 


If  bhe  wavelength  and  time-dependent  source  term  S(X,t) 

and  F2(T-t)  is  the  atmospheric  transmission  term  I(X,T-t), 
then  the  convolution  process  becomes 


M(X,t)  = 


S(X,t)  I(X,T-t)dt 


which  for  this  comparison  results  in 


a a r -a  x -a,x1 


Assuming  a detector  response  R(X)  such  that 


R(X)  = 1:  1 X 1 >^2 


and  R(X)  = 0:  elsewhere, 

integration  over  wavelength  yields  the  detected  signal 


f,  - 

• 

4 

-a,T  -a,T 

X"  - X, 

1 1 
e -e 

2 1 

For  check-out  purposes,  Xj^  = .4278  microns,  X2  = 1.07  microns,  thus 


D(X)  = 5.322X10 

.-5 


-5 


■ -io\  -io\' 

e -e 


which  has  a maximum  of  3.71X10  at  25.6  ysec.  The  results  of  the  analytic 
solution  are  shown  as  a continuous  curve  in  Fig.  20.  The  discrete  results 
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Fig.  20.  Comparison  of  TMTAU,  POLYN,  and  TMTRFOl  Calculations  for 
Check-Out  Problem  Using  an  Analytic  Source  Term  and  an 
Analytic  Transmission  Function 
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of  TMTAU  and  TMTRFOl  are  plotted  at  the  evaluated  times;  POLYN  results  are 
plotted  at  the  logarithmic  mean  times  on  the  same  graph.  No  direct  radiation 
contribution  was  assumed  in  the  calculations.  Results  from  each  of  the  three 
computer  procedures  are  shown  only  where  the  discrepancy  between  the  analytic 
so2,ution  and  the  results  from  each  procedure  were  noticeable.  Apparent 
errors  in  the  results  from  TMTAU,  TMTRFOl,  and  POLYN  for  early  times  can  be 
explained  by  the  differences  in  input  data  for  those  codes.  TMThU  was  run 


TMTRFOl  and  POLYN  were 
.-3 


-7  -3 

with  the  source  defined  for  times  10  < t < 10 

run  with  the  source  defined  for  times  10  ^ < t < 10  Since  the  convolution 
integral  as  evaluated  above  had  a lower  time  limit  of  zero,  any  other  lower 
bound  will  introduce  an  error.  That  error  is  apparent  when  the  integral 


M’(A,t)  e 


S(X,t)  I (X,T-t)dt 


0 

is  evaluated  and  compared  to  the  convolution  solution  for  a lower  time  bound 
of  0.0.  The  results  are: 


M(X,t)  = 


^0^0 

(a^-ap 


-aiT 

e -e 


-C(  t] 


M’a.x) 


(ai-a^) 


r-v 


-»T 

- e e 


where  the  vanishing  error  is  e . The  error  is  vanishing  since 

-V  "V„  -V  -V  S-'i^^o  -"i"  “i 

e -e  - e -e  e = e for  — = 10, 

— 8 — 3 

Tq  < T and  a^T  >5.  A run  of  TMTRFOl  for  times  10  ^ ^ 10  resulted  in  ... 

absolute  error  of  less  than  .05%  at  one  microsecond.  A comparison  of  the 
POLYN  and  TMfRFOl  results  for  times  10  ^ < t < 10  ^ and  the  analytic  formula 


M'(X,t)  = 


S(X,t)  I(X,T-t)dt, 
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where  = 10  is  shown  In  Fig.  21.  The  variance  between  the  TMTAU, 
TMTRFOl  and  POLYN  results  and  the  analytic  results  is  only  aesthetically 
unpleaslng  since  it  will  only  be  apparent  in  sources  which  have  signi- 
ficant values  at  zero  source  time.  This,  of  course,  does  not  occur  in 

realistic  sources  and,  as  evidenced  by  the  TMTRFOl  results  where  the 
—8  —3 

times  were  10  ^ t ^ 10  5 ihe  procedure  results  can  adequately  define 

the  convolution  process  if  the  parametric  data  adequately  defines  the 
source  and  transmission. 

One  of  the  basic  uncertainties  in  the  use  of  TMTAU  has  been  in 
the  results  of  the  integration  over  wavelength  at  a given  time  which  was 
based  on  the  use  of  source  and  atmospheric  transport  data  for  only  three 
wavelengths,  0.4278  pm,  0.75  pm,  and  1.07  um.  This  uncertainty  resulted 
in  the  development  of  the  TMTRFOl  procedure  which  uses  curve-fit  coeffi- 
ci i ,.ts  data  for  more  than  three  wavelengths  in  the  convolution  integral. 
The  uncertainty  in  the  TMTAU  data  has  also  resulted  in  a renewed  interest 
in  the  POLYN  procedure  which  also  has  the  capability  of  doing  the  convo- 
lution Integral  for  more  than  three  wavelengths. 

The  TMTRFOl  procedure  differs  from  TMTAU  only  in  that  it  provides 
for  the  use  of  more  than  three  wavelengths.  The  wavelength  dependence 
of  the  power  at  the  detector  at  a given  time  is  expressed  in  terms  of 
the  function 

(a  (t)  + a (t)  X + a (t)  X^  • • + a X'^) 

P(t,X)  = e ^ 2 3 

in  TMTRFOl  where  n is  the  number  of  wavelengths  for  which  the  source  and 
atmospheric  transmission  data  are  defined.  The  integral  over  wavelength 
is  then  evaluated  numerically  with  use  of  the  equation 

28 

P(t)  = I P(t,X  ) R(X  )AX 
1>1 

in  both  codes  with  an  Integration  Interval  of  AX  = (X^-Xj^)/28  where  X^ 
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Comparison  of  POLYN  and  TMTRFOl  Results  for  Test  Problem,  Lower  Time  Value 
for  Source  Function  • 10“6  sec  in  all  Solutions 


r 


and  are  the  upper  and  lower  bounds  of  the  wavelength  range  = 

1.1  pm  and  X.  = 0.40  pm)  and  the  X are  the  midpoints  of  each  of  the 

Xr  1 

28  subintervals  in  wavelength  into  which  the  interval  is  divided. 

In  POLYN  the  integration  is  performed  over  wavelength  using  AX  values 
input  for  each  wavelength  for  which  the  receiver  response  data  are  input. 

A test  problem  was  run  the  POLYN  and  TMTRFOl  codes  using  trans- 
mission data  obtained  by  interpolation  for  14  wavelengths  from  the  curve- 
fit  coefficient  data  for  the  40-km  meteorological  range  atmosphere  at 
X = 0.4278  pm,  0.75  pm,  and  1.07  pm.  The  RADFLO  source  term  for  a 2-KT 
source  at  0.1-km  altitude  with  the  silicon  response  function  (Fig.  9) 
for  the  detector  response  was  also  used.  The  wavelength  dependence  of 
the  power  at  the  detector  at  satellite  altitude  (35800  km)  for  a look 
angle  of  0°  as  computed  by  each  of  the  convolution  codes  is  given  in 
Tables  XVII,  XVIII,  and  XIX  for  retarded  times  of  128.6  psec,  5.984  msec, 
and  52.64  msec,  respectively.  The  data  given  in  Tables  XVII,  XVIII,  and 
XIX  for  the  TMTAU  procedure  were  obtained  from  evaluation  of  the  equation 

(a  (t)  + a (t)X  + a (t)X^) 

P(t)  = e ^ 

where  the  a^  coefficients  were  obtained  from  the  TMTAU  output  for  the 
indicated  times.  The  a^  coefficients  were  fitted  to  the  TMTAU  results 
for  wavelengths  of  0.4278,  0.75,  and  1.07  p.  Also  shown  in  the  tables 
is  the  Integral  over  wavelengths  for  each  of  these  times.  It  is  noted 
that  POLYN  and  TMTRFOl  produce  results  as  a function  of  wavelength  that 
differ  mainly  only  in  the  second  digit  at  the  times  shown.  It  is  noted 
that  the  wavelength  dependence  of  the  power  as  given  by  TMT.AU  differs 
significantly  from  that  given  by  POT.YN  and  TMTRFOl.  The  times  useJ  • 

Tables  XVII,  XVIII,  and  XIX  represent  times  near  first  maximum,  mi 

and  second  maximum  of  the  power  history  curve.  The  ratios  i t P'  . 

TMTRFOl  and  TMTAU  to  TMTRFOl  for  the  Integral  over  waveltnxt' 
times  are  shown  in  Table  XX. 
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TABLE  XVII.  CONVOLUTION  INTEGRAL  AS  A FUNCTION  OF  WAVELENGTH 
AT  A TIME  OF  128.6  psec 

2 

POWER  (watts/cm  -p) 


Wavelength 

(microns) 

TMTRFOl 

POLYN 

TMTAU 

.4278 

1.847+2 

1.860+-2 

1.793+2 

.4675 

4.156+2 

4.249+2 

2.400+2 

.5052 

5.497+2 

5.660f2 

2.998+2 

.5476 

3 , 792+2 

3.915+2 

3.617+2 

.5983 

4.525+2 

4.654+2 

4.144+2 

.6539 

3.213f2 

3.278+2 

4.311+2 

.7085 

3.60^2 

3.648f2 

4.007+2 

.7693 

2.239+2 

2.243+2 

3.241+2 

.8364 

2.463+2 

2.453+2 

2.187+2 

.9004 

9.058+1 

8.988+1 

1.285f2 

.9613 

9.112+1 

9.066+1 

6.728+1 

1.0059 

6.638+1 

6.626+1 

3.837+1 

1.0400 

3.803+1 

3.816+1 

2.37&fl 

1.0700 

1.555+1 

1.568+1 

1.504+1 

2 

Integral  (watts/cm  ) 

1.651+2 

1.752+2 

1.70f2 

TABLE  XVIII.  CONVOLUTION  INTEGRAL  AS  A FUNCTION  OF  WAVELENGTH 
AT  A TIME  OF  5.984  msec 


2 

POWER  (watts/cm  -y) 


Wavelength 

(microns) 

TMTRFOl 

POLYN 

TMTAU 

.4278 

3.715+0 

3.812+0 

3.572+0 

.4675 

8.309+0 

8.663+0 

5.760+0 

.5052 

1.093+1 

1.148+1 

8.614+0 

.5476 

1.950H 

1.978+1 

1.277+1 

.5983 

2.312+1 

2.336+1 

1.881+1 

.6539 

3.234+1 

3.172+1 

2.595+1 

.7085 

3.615+1 

3.502+1 

3.203+1 

.7693 

4.745+1 

4.469+1 

3.582+1 

.8364 

5.204+1 

4.872+1 

3.487+1 

.9004 

2.823+1 

2.618+1 

2.935+1 

.9613 

2.836+1 

2.634+1 

2.181+1 

1.0059 

4.050+1 

3.760+1 

1.615+1 

1.0400 

2.318+1 

2.162+1 

1.225+1 

1.0700 

9.466+0 

8.874+0 

9.291+0 

Integral  (watts 

/cm^) 

1.941+1 

1.855+1 

1.51+1 

r 


i 

I 


f 

[ 

TABLE  XIX.  CONVOLUTION  INTEGRAL  AS  A FUNCTION  OF  WAVELENGTH 
AT  A TIME  OF  52.64  msec 

2 

POWER  (watts/cm  -p) 


Wavelength 


(microns) 

TMTRFOl 

POLYN 

TMTAU 

.4278 

3.368+2 

3.315+2 

3.279+2 

.4675 

7.534+2 

7.535+2 

4.433+2 

.5052 

9.912+2 

9.990f2 

5.615+2 

.5476 

8.032+2 

8.137+2 

6.910+2 

.5983 

9.524+2 

9.609+2 

8.171+2 

.6539 

7.164+2 

7.155+2 

8.875+2 

.7085 

8.007+2 

7.899+2 

8.685+2 

.7693 

4.882+2 

4.760+2 

7.519+2 

.8364 

5.352+2 

5.189+2 

5.538+2 

.9004 

2.402+2 

2.323+2 

3.584+2 

.9613 

2.412+2 

2.337+2 

2.080+2 

1.0059 

2.414+2 

2.358+2 

1.288+2 

1.0400 

1.382+2 

1.350+2 

8.533+1 

1.0700 

5.641+1 

5.565+1 

5.746+1 

2 

Integral  (watts/cm  ) 

3.540+-2 

3.656+2 

3.62+2 

TABLE  XX.  RATIO  OF  THE  INTEGRAL  OVER  WAVELENGTH  AS  GIVEN  BY 
THE  POLYN  AND  TMTAU  CODES  TO  THAT  GIVEN  BY  TMTRFOl 
CODE  FOR  TIMES  OF  128.6  psec,  5.984  msec,  and  52.64 
msec 
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A comparison  of  the  time  histories  of  the  power  at  the  detector 
as  given  by  the  POLYN,  TMTRFOl,  TMTAU  codes  Is  given  In  Table  XXI.  Con- 
sidering the  differences  In  the  methods  used  In  each  code,  there  Is 
reasonably  good  agreement  between  the  results  of  the  three  codes.  (Note 
that  the  powers  as  given  by  POLYN  and  TMTRFOl  have  to  be  divided  by  a 
factor  of  10^^  to  put  them  In  the  same  units  as  that  given  by  TMTAU.) 

In  light  of  the  results  given  In  Tables  XVII  through  XXI,  It  Is  suggested 
that  POLYN  or  TMIRFOl  be  used  In  the  future  Instead  of  TMTAU  to  avoid 
the  uncertainties  In  the  wavelength  Integral  for  each  time  when  only 
three  wavelengths  are  used  to  represent  both  the  source  and  atmosphere 
transmission  data. 

6.3  Conversion  of  TMTAU  Calculations  for  Synchronous  Satellite 
Altitude  to  Other  Satellite  Altitudes 

POLO  calculations  were  run  to  determine  If  the  time-dependent 
scattered  fluxes  at  receiver  altitudes  of  20186.8  km  and  93621.626  km 
could  be  determined  from  the  calculations  run  for  a detector  altitude 
of  35800  km.  The  POLO  calculations  were  run  for  a 0.75  vi  wavelength 
point  Isotropic  source  positioned  at  a height  of  0.5  km  above  the 
ground.  The  ground-level  meteorological  range  was  25  km.  Receivers 
were  positioned  at  altitudes  of  20186.8  km,  35800  km,  and  93621.626 
km.  Scattered  fluxes  as  a function  of  time  after  arrival  of  the  direct 

_3 

flux  were  computed  for  50  time  Intervals  between  0 and  1.0X10  seconds. 
The  detectors  were  posltlon?.d  for  each  detector  altitude  at  polar  angles, 
, between  a radial  through  the  source  and  the  line  connecting  the 
source  and  receiver,  of  0,  30,  60,  and  90  degrees.  The  POLO  geometry 
Is  shown  In  Fig.  22.  R^,  R^,  and  R^  are  the  distances  from  the  source 

to  the  satellite  altitudes  of  20186.8  km,  35800  km,  and  93621.626  km, 
respectively.  In  the  direction  denoted  by  6 . The  angle  between  the 

MU 

source-receiver  axis  and  the  earth  radial  through  the  satellite  Is 
denoted  as  8.  The  uncolllded  flux  at  each  receiver  for  a unit  point 
Isotropic  source  Is  given  by  the  expression. 
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TABLE  XXI.  COMPARISON  OF  POLYN,  TMTRFOl  AND  TMTAU 
CALCULATIONS 


F 

L. 


POWER  (watts/cm  ) 
POLYN*  TMTRFOl* 


IMTAU 


3.260E-06 

6t392E-06 

1 

OOOE 

319E 

00 

00 

1.173E  00 

1 .14E-10 

7ta70E-06 

1 

320E 

00 

1*230E  00 

1 . 21E-10 

1 

275E 

00 

1«191E  00 

1 .18F-10 

l«baSE*05 

1 

562E 

00 

1»860E  00 

1 .9t)E-lO 

2«l<»0E-05 

2 

236E 

00 

2«440E  00 

2 .53E-10 

3i006E-05 

4 

495E 

00 

A»858r  00 

6. 7flE-in 

3*289E**05 

7 

421E 

00 

7t463E  00 

7.49E-10 

3.7<»lE-05 

1 

002E 

01 

1^27AE__Q1... 

l.lie-09 

<».30hE-05 

2 

091E 

01 

3«268E  01 

2.96E-09 

<n7‘»8E“05 

3 

909E 

01 

h>72hE  01 

4 .87E-09 

<»*930E-05 

5 

S22E 

01 

7>066E  01 

7 . B2E-n9 

5»‘»20E-05 

8 

588E 

01 

9t846E  01 

1 . llE-08 

StS'tlE-OS 

1 

191E 

02 

li3<»lE  02 

1.44E-08 

7.076E-05 

1 

601E 

02 

_l!74.8E.  02 

1 .82E-0B 

7*69bE-05 

1 

S36£ 

02 

lt777E  02 

1 .84E-08 

Bt366E-05 

1 

8&8E 

02 

lt784E  02  - 

1 .84E-08 

9*0ii2E*‘05 

1 

8&8E 

02 

1*769E  02 

1.83E-08 

1 .lOlE-04 

1 

818E 

02 

1.700E  02 

1 . 76E-08 

1 

776E 

02 

lt683E  02 

1 . 74E-08 

lt286E-0^ 

1 

752E 

02 

1.651E  02 

1 . 70E-08 

l*559E-0-» 

1 

68SE 

02 

1<557E  02 

1 .60E-08 

lt8‘t7E“04 

1 

581E 

02 

li463E  02 

1 .51E-08 

2«l‘»lE-0% 

1 

4S3E 

02 

li303E  02 

1 .34E-08 

2><»08E«04 

1 

340E 

02 

1.223E  02 

1 .26E-08 

2>676E-0^ 

1 

2S0E 

02 

lil44E  02 

1.17E-08 

2«939E*'0‘r 

1 

170E 

02 

K078E  02 

1 . 11 E-06 

3»206E»04 

1 

127E 

02 

1*0S2E  02 

1 . 08F-C8 

3»<»73E-0^ 

1 

106E 

02 

1«03SE  02 

1 . 06F-08 

3i7<rOE-0*» 

1 

032E 

02 

9*200E  01 

9 .41E-09 

4«007E-04 

9 

70<»E 

01 

9t217E  01 

9 .43E-09 

%«407£"0<» 

9 

186E 

01 

8>13SE  01 

8.31E-09 

4«67<tE*0H 

8 

362E 

01 

7*693E  01 

7 .35E-09 

5*3*»aE-0-» 

7 

667E 

01 

6>8S1E  01 

6 . 98E-09 

6«008E-0^ 

6 

4B6E 

01 

__S*497E  01 

5.58E-09 

6>675E-0<t 

5 

6t>7E 

01 

5«310E  01 

5 .38P-P9 

7»342E-04 

5 

267E 

01 

4*709E  01 

4 .76E-09 

B*008E*0% 

4 

8i»4E 

01 

4*S44E  01 

4.58E-09 

S».3«»1E-04» 

** 

5S0E 

0l 

4*147E  01 

4 . 16E-09 

1.067E-03 

4 

331E 

0l 

4*139E  01 

4 . 14E-P9 

1 •201E‘03 

4 

381E 

01 

4*25%E  01 

4 .25E-09 

1*33<»E*03 

4 

314E 

01 

4i017E  01 

4 .OlE-09 

1 (bOOE-OS 

4 

035E 

01 

3>737E  01 

3.71E-09 

li866E-03 

3 

695E 

01 

3«363E  01 

3.31F-09 

2*133E-03 

3 

278E 

01 

2«959E  01 

? . 38E-09 

2*399E-03 

2 

870E 

01 

2>S96E  01 

2 . 48E-P9 

2*66bE*03 

2 

509E 

01 

2*278E  01 

2.13E-09 

2*931E-03 

2 

2S9E 

01 

2*121E  01 

1 .9'5F-09 

3<197E-03 

2 

130E 

Ol 

2*03bE  01 

1 .84E-09 

3«462E-03 

2 

031E 

01 

1i941E  01 

1 . 73E-U9 

3>72llE-03 

1 

954E 

01 

1>896E  01 

1 . 66E-09 

* POLYN  and  TMTRFOl  results  are  10 
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POWER 


i 


TABLE  XXI.  (Continued) 


POWER  (watts/cm^) 


TIME 

SECONDS 

POLYN* 

TMTRFOl* 

TMTAU 

3 

• 994E'-03 

1 

.897E 

01 

1 

• 839E 

01 

1.59E-C9 

•392E-03 

1 

• a49E 

01 

i 

• Sl4£ 

01 

1 . 54E-09 

•69BE*03 

1 

• a33E 

01 

1 

• 816E 

01 

1.52E-09 

5 

•321E-03 

1 

• a02E 

01 

i' 

.774E 

01 

1 . 41E-09 

5 

•9B4E-03 

1 

• 895E 

01 

1 

• 941E 

01 

1 .51E-09 

6 

•647E*03 

2 

• 06BE 

01 

2 

.207E 

01 

1 . 71E-09 

7 

•310E-03 

2 

• 288E 

01 

2 

• 376E 

01 

t .66E-09 

7 

.972E-03 

2 

.791E 

01 

3 

• 240E 

01 

2.69E-09 

9 

•296E-03 

3 

• 619E 

01 

3 

• 943E 

01 

3.27E-09 

1 

•062E-02 

4 

• 627E 

01 

5 

• 252E 

01 

4 . 61E-09 

1 

•194E-02 

6 

• 187E 

01 

6 

• 951E 

01 

6 . 39E-09 

1 

•326E-02 

8 

• 086E 

01 

8 

*866E 

01 

8 . 46E-09 

1 

•590E-02 

1 

• 066E 

02 

1 

• 184E 

02 

1.15E-08 

1 

• il5<»E-02 

1 

• 426E 

02 

1 

.570E 

02 

1 .56E-08 

2 

•llBE-02 

1 

.815E 

02 

1 

• 912E 

02 

1.93E-08 

2 

•3O1E-02 

2 

• 152E 

02 

2 

• 209E 

02 

2.23E-06 

2 

• 644E*02 

2 

• 436E 

02 

2 

t451E 

02 

2.50E-08 

3 

•169E-02 

2 

• 767E 

02 

2 

•832E 

02 

2 .B9E-08 

3 

• <»32E-02 

3 

• C36E 

02 

2 

*999E 

02 

3. 07E-0B 

3 

•694E-02 

3 

• 225E 

02 

3 

• 153E 

02 

3.23E-08 

3 

•996E-02 

3 

• 364E 

02 

3 

• 261E 

02 

3. 34E-06 

<» 

• 3<»9E"02 

3 

• 485E 

02 

3 

• 381E 

02 

3 . 46E-ns 

<» 

•611E-02 

3 

.572E 

02 

3 

• 427E 

02 

3 .51E-08 

5 

•264E»02 

3 

• 6ti6E 

02 

3 

• S40E 

02 

3 .62E-08 

9 

•917E-02 

3 

• 745E 

02 

3 

• 598E 

02 

3.68E-08 

6 

•369E-02 

3 

.7592 

02 

3 

• SbiE 

02 

3 .64E-08 

7 

>220E-02 

3 

• 6S1E 

02 

3 

.41<»E 

02 

3.48E-08 

7 

.H70E-02 

3 

• 424E 

02 

3 

• 142E 

02 

3.20E-08 

9 

. 168E-02 

2 

• 849E 

02 

2 

• 412E 

02 

2.35E-08 

1 

•046E-01 

2 

• 130E 

02 

1 

• 799E 

02 

_ 1.60E-0B 

1 

•17bE-01 

1 

• 63BE 

02 

-A 

•>69E. 

02 

1.15E-08 

1 

.305E-01 

1 

• 346E 

02 

1 

• 244E 

02 

9.09E-09 

1 

•SbbE-Ol 

1 

• 058E 

02 

9 

• 146E 

01 

6 .35E-09 

1 

• b9oE'*01 

8 

.4*»0E 

01 

»^915£_Q1 

5.51E-09 

1 

f B27E-01 

7 

t364E 

01 

6 

• 972I_ 

01 

4 .87E-09 

1 

•937E-01 

6 

• 495E 

01 

-i 

«1S9E 

Pi 

4 . 3bE-09 

2 

•218E-01 

S 

• 49SE 

01 

4 

«960E  01 

3.6nE-09 

2 

•6l0E*0l 

4 

• 3S6E 

01 

_i 

• 786E 

_Pi__ 

2.91E-09 

2 

•870e*01 

3 

• 579E 

01 

3 

• 36SE 

01 

2.58E-09 

3 

•131E-01 

3 

.U9E 

01 

2 

• 878E 

01 

2.31E-09 

3 

•392E-01 

2 

. 775E 

01 

2 

•597E 

01 

2.13E-09 

3 

•653E-01 

2 

• 491E 

01 

2 

• 316E 

01 

1.95E-09 

3 

• 9l*»E-01 

2 

• 266E 

01 

2 

• 14AE 

01 

1 .83F-09 

4 

•30bE-01 

2 

• 031E 

01 

1 

.924E 

01 

1 . 66E-09 

*» 

I697E-01 

1 

• b12E 

01 

1 

t670E 

01 

1 . 50E-09 

5 

•219E-01 

1 

• 605E 

01 

1 

• 486E 

01 

1 .36E-09 

5 

•871E-01 

1 

.427E 

01 

A 

• 318E 

01 

1 . 24E-09 

*POLYN  and  TMTRFOl  results  are  10^®  X POWER 
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Fig.  22.  POLO  Geometry 


Ijj(R)  - e'^cose/ATTR^ 


where  p is  the  optical  distance  in  the  atmosphere  along  the  path  of 
the  direct  radiation.  Therefore,  the  uncolllded  flux  at  R^  in  terms 
of  the  uncolllded  flux  at  Rg  is  given  by 


W 


R?  cose 

^ coieT 


This  results  from  the  fact  that  the  atmospheric  attenuation  along  the 
path  of  the  direct  radiation  in  the  atmosphere  is  the  same  no  matter 
where  the  receiver  is  positioned  above  the  atmosphere.  The  angles  By  > 

Li 

Bg,  By  are  given  by 


° ®SL~®L’  ^ “ ®Sl“®V 


Table  XXll  lists  the  value  of  0 , 6_,  and  0 , and  B,  , B«»  B,,  for 

L D V L S V 

the  three  satellite  altitudes  when  0g^^  «=  0,  30,  60,  and  90  degrees. 

From  the  geometry,  one  can  obtain  the  values  of  the  quantities 
2 2 2 2 

RgCOsBj^/Rj^cosBg  and  RgCosBy/RyCosBg.  These  values  are  listed  in  Table 

XXIII  as  a function  of  the  look  angle  0 . In  each  of  the  50  retarded 

SL 

time  intervals,  the  ratios  of  the  scattered  fluxes  that  were  computed 
for  the  satellite  altitude  of  20186.8  km  to  that  computed  for  the 
satellite  altitude  of  35800  km  were  found  to  agree  with  the  values  of 
the  factors  listed  in  Table  XXIII  for  the  first  three  significant  digits 
for  the  four  different  look  angles.  Similarly,  the  ratios  of  the  scat- 
tered fluxes  as  a function  of  retarded  time  for  the  satellite  altitude 
of  93621.626  to  that  for  the  satellite  altitude  of  35800  km  were  found 
to  also  agree  with  the  values  listed  in  Table  XXIII  within  the  first 
three  significant  figures.  In  fact,  nearly  all  of  the  ratios  computed 
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TABLE  XXII.  VALUES  OF  0^,  6^,  0^,  0^,  0^,  AND  0^  FOR  EACH  SATELLITE 
ALTITUDE  AS  A FUNCTION  OF  0 

SL 


®SL 

(deg) 

Altitude= 

20186.8  km 

Altitude= 

93621.626  km 

Altitude- 

35800  km 

®L 

(deg) 

^L 

(deg) 

®V 

(deg) 

®V 

(deg) 

®S 

(deg) 

^S 

(deg) 

0 

0 

0 

0 

0 

0 

0 

30 

23.10437 

6.89563 

28.17228 

1.82772 

25.66327 

4.33673 

60 

47.99766 

12.0023 

56.83321 

3.16679 

52.47412 

7.52588 

90 

76.10626 

13.8937 

86.34269 

3.65731 

81.30145 

8.69855 

TABLE  XXIII.  FACTORS  FOR  CONVERTING  DIRECT  FLUXES  FROM  SYNCHRONOUS 
SATELLITE  ALTITUDE  TO  ALTITUDES  OF  20186.8  KM  AND 
93621.626  KM 


®SL 

Rgcos0j^/R^cos0g 

Altitude-20186.8  km 

2 2 
RgCOS0y/R^COS0g 

Altitude-93621. 626  km 

0 

3.145 

0.1462 

30 

3.057 

0.1500 

60 

2.833 

0.1608 

90 

2.567 

0.1762 
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for  the  scattered  fluxes  were  found  to  differ  by  no  more  than  0.005  for 
the  satellite  altitude  of  20186.8  km  and  no  more  than  0.0005  for  the 
satellite  altitude  of  93621.626  km  from  the  values  listed  in  Table  XXIII. 

From  the  results  of  the  POLO  calculations,  it  is  concluded  that 
direct  and  time-dependent  scattered  fluxes  for  satellite  altitudes  of 
20186.8  km  and  93621.626  km  can  be  obtained  with  an  error  of  less  than 
1.0%  from  the  data  previously  generated  with  POLO  for  a satellite  alti- 
tude of  35800  km  (Ref.  14).  These  data  are  stored  on  Tape  No.  2696 
(labeled  POLODCBINE ,DSN  = POLOD)  in  the  order  shown  In  Table  XXIV.  The 
cumulative  distributions  of  the  scattered  fluxes  were  curve-fitted  with 
use  of  the  expression. 


i<6 


Is(t)dt 


^1  - ? 


-b  .T 

a^e  ^ + 


j<3 

-I 


-b.(T-To) 


where  - maximum  cumulated  scattered  flux  from  photons  that  have 
not  undergone  ground  reflection, 

K2  “ maximum  cumulated  scattered  flux  from  photons  that  have 

undergone  both  ground  reflection  and  air  scattering  minus 

Tq  “ retarded  time  at  which  ground-reflected  photons  begin  to 
arrive  at  the  receiver. 


and  a^,  b^,  a^ , and  b^  are  constants  obtained  by  curve-fitting  the  cumu- 
lative flux  values.  The  curve-fit  coefficients  for  the  synchronous 
satellite  altitude  (35800  km)  are  stored  on  Tape  No.  801  (labeled 
POLODCVFITSC,DSN  “ POLOFT)  according  to  the  sequence  shown  in  Table 
XXIV.  The  curve-fit  data  on  Tape  No.  801  are  for  look  angles  of  0, 

0,  30,  45,  60,  70,  75,  80,  and  90°. 


LS 


To  obtain  a power  history  at  satellite  altitudes  of  20186.8  km 
or  93621.626  km  for  a given  source  power  history,  it  is  suggested  that 
the  POLO  data  on  Tape  No.  801  be  used  in  the  TMTAU  procedure.  The 
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TABLE 

XXIV.  SEQUENCE 
NO.  2656 

OF  POLO  1 
(I ABELED 

OUTPUT  DATA  SETS 
POLODCBINE,  DSN 

STORED  ON  TAPE 
- POLOD) 

Met. 

Source 

Source 

Sequence 

Problem 

Range 

Alt. 

Wavelength 

No. 

No. 

(km) 

(km) 

(ym) 

1 

43031 

3 

0.001 

0.4278 

2 

75031 

3 

0.001 

0.75 

3 

107031 

3 

0.001 

1.07 

4 

43032 

3 

0.5 

0.4278 

5 

75032 

3 

0.5 

0.75 

6 

107032 

3 

0.5 

1.07 

7 

43033 

3 

3.0 

0.4278 

8 

75033 

3 

3.0 

0.75 

9 

107033 

3 

3.0 

1.07 

10 

43034 

3 

10.0 

0.4278 

11 

75034 

3 

10.0 

0.75 

12 

107034 

3 

10.0 

1.07 

13 

43101 

10 

0.001 

0.4278 

14 

75101 

10 

0.001 

0.75 

15 

107101 

10 

0.001 

1.07 

16 

43102 

10 

0.5 

0.4278 

17 

75102 

10 

0.5 

0.75 

18 

107102 

10 

0.5 

1.07 

19 

43103 

10 

3.0 

0.4278 

20 

75103 

10 

3.0 

0.75 

21 

107103 

10 

3.0 

1.07 

22 

43104 

10 

10.0 

0.4278 

23 

75104 

10 

10.0 

0.75 

24 

107104 

10 

10.0 

1.07 

25 

43251 

25 

0.001 

0.4278 

26 

75251 

25 

0.001 

0.75 

27 

107251 

25 

0.001 

1.07 

28 

43252 

25 

0.5 

0.4278 

29 

75252 

25 

0.5 

0.75 

30 

107252 

25 

0.5 

1.07 

31 

43253 

25 

3.0 

0.4278 

32 

75253 

25 

3.0 

0.75 

33 

107253 

25 

3.0 

1.07 

34 

43254 

25 

10.0 

0.4278 

35 

75254 

25 

10.0 

0.75 

36 

107254 

25 

10.0 

1.07 
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results  so  obtained  from  TMTAU  for  a satellite  altitude  of  35800  km  can 
then  be  converted  to  either  the  20186.8  km  or  93621.626  km  satellite 
altitudes  with  the  use  of  the  correction  factors  listed  In  Table  XXV. 

It  must  be  pointed  out  that  only  the  angle  0 remains  a constant  when 
converting  data  for  a satellite  altitude  of  35800  km  to  other  satellite 
altitudes;  that  Is,  the  look  angles  6 ^ and  6 are  different  for  each 

IjV 

source  altitude  and  look  angle  0 used  In  the  POLO  calculations. 

L<d 

An  examination  of  Table  XXV  reveals  that  the  correction  factors 
and  look  angles  for  each  satellite  altitude  (20186  , and  93621.626  km) 
do  not  vary  significantly  with  source  altitude.  It  app-  ars  for  all 
source  altitudes  that  one  can  assume  that  Qfc  “ 

JLj^  UU  IrfV 

correction  factors  as  a function  of  0 for  a source  altitude  of  3.0  km 
can  be  used  for  all  source  altitudes  without  Introducing  any  significant 
error  In  the  results  obtained. 
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TABLE  XX'/ 


FACTORS  FOR  CONVERTING  TMTAU  AND  POLO  CALCULATIONS  FOR  A 
SATELLITE  ALTITUDE  OF  35800  KM  TO  SATELLITE  ALTITUDES  OF 
20186.8  KM  AND  93621.626  KM 


Source 

Height 

(km) 

®LS 

35800km  to 

20186.8km 

35800km  to  93621.626km 

6LL(deg) 

FACTOR 

FACTOR 

0.001 

0 

0.00000 

3.14520 

0.00000 

0.146224 

30 

29.99999 

3.05705 

30.00000 

0.150050 

45 

44.99999 

2.95672 

45.00000 

0.154675 

60 

59.99999 

2.83285 

60.00000 

0.160832 

70 

69.99999 

2.74374 

70.00000 

0.165619 

75 

74.99999 

2.69870 

75.00000 

0.168169 

80 

79.99999 

2.65396 

80.00000 

0.170798 

90 

89.99999 

2.56699 

90.00000 

0.176214 

0.5 

0 

0.00000 

3.14527 

0.00000 

0.146216 

30 

29.99971 

3.05711 

30.00024 

0.150048 

45 

44.99960 

2.95676 

45.00033 

0.154674 

60 

59.99955 

2.83288 

60.00038 

0.160831 

70 

69.99954 

2.74376 

70.00040 

0.165618 

75 

74.99954 

2.69872 

75.00041 

0.168168 

80 

79.99955 

2.65397 

80.00041 

0.170798 

90 

89.99820 

2.56702 

90.00097 

0.176215 

3.0 

0 

0.00000 

3.14560 

0.00000 

0.146209 

30 

29.99823 

3.05738 

30.00144 

0.150037 

45 

44.99764 

2.95698 

45.00197 

0.154665 

60 

59.99732 

2.83302 

60.00231 

0.160825 

70 

69.99727 

2.74385 

70.00243 

0.165614 

75 

74.99728 

2.69879 

75.00245 

0.168166 

80 

79.99732 

2.65402 

80.00245 

0.170796 

90 

89.98921 

2.65672 

90.00584 

0.176218 

10.0 

0 

0.00000 

3.14656 

0.00000 

0.146174 

30 

29.99410 

3.05817 

30.00481 

0.150007 

45 

44.99213 

2.95758 

45.00658 

0.154641 

60 

59.99108 

2.83342 

60.00771 

0.160809 

70 

69.99089 

2.74411 

70.00809 

0.165604 

75 

74.99093 

2.6989L 

75.00817 

0.168158 

80 

79.99106 

2.65415 

80.00818 

0.170792 

90 

89.96404 

2.56764 

90.01948 

0.176227 
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VII.  THERMAL  RADIATION  TRANSPORT  IN  CLOUDY  ATMOSPHERES 


A series  of  POLO  Monte  Carlo  calculations  were  run  during  the 
reporting  period  to  study  the  effects  of  cloud  layers  on  the  time- 
dependent  trans'  . .t  of  the  radiation  emitted  from  point  isotropic 
sources  to  detecto/s  positioned  at  synchronous  satellite  altitude. 

The  model  atmosphere  used  in  the  Monte  Carlo  calculations  had 
a ground-level  meteorological  range  of  10  km.  The  size  distribution 
used  to  characterize  the  aerosols  was  the  Haze  C distribution  (Ref.  17) 
for  continental  haze.  POLO  problems  were  run  for  two  different  cloud 
models.  In  the  first  POLO  problems,  a stratus  cloud  with  an  optical 
thickness  of  approximately  14  mean-free-paths  at  all  wavelengths  was 
positioned  between  0.5  km  and  1.0  km  altitude.  In  the  next  set  of 
POLO  problems,  an  altostratus  cloud  layer  with  approximately  a ?0 
mean-f ree-path  thickness  at  all  wavelengths  was  positioned  between 
3 km  and  4 km  altitude. 

The  particle  size  distribution  N(r)  used  to  describe  the  cloud 
models  is  given  by  the  equation  (Ref.  21) 


N(r)  = A r e 


•Br' 


The  total  number  Nq  of  particles  per  cm  is  given  by  integrating  the 


size  distribution  over  all  radii. 


If  one  assumes  that  r . =0  and 

min 


r = then 
max 


YB 


(ot+l)/Y 


N(r)dr 


where  F(x)  is  the  gamma  function.  If  N^,  a,  B,  and  y are  known,  then 
A can  be  determined  from  the  above  equation.  For  the  stratus  cloud: 
a = 3,  B = 0.666,  y “ 1,  » 250,  and  the  size  range  was  0.4  to  30  ym 

For  the  altostratus  cloud:  a - 5,  B - 1.111,  y - 1,  N^  - 400,  and  the 
size  range  was  0.2  to  18  ym.  The  droplets  for  each  model  were  assumed 
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to  be  water.  The  index  of  refraction  for  water  at  the  wavelengths 
used  In  the  POLO  calculations  was  taken  from  Ref.  22. 


Point  isotropic  sources  were  positioned  at  altitudes  of  0.1, 
0.5,  3,  and  10  km.  The  POLO  problems  were  run  for  source  wavelengths 
of  0.4278,  0.75,  and  1.07  microns.  Detectors  were  positioned  at 
satellite  look  angles  of  0,  45,  75,  and  90  degrees.  The  ground  surface 
was  assumed  to  be  a Lambert  reflector.  The  cloud  layer  was  taken  to  be 
a scattering  medium  and  not  as  a reflection  surface.  The  POLO  calcula- 
tions were  run  for  a satellite  altitude  of  35800  km. 

The  POLO-computed  direct  and  time-dependent  scattered  fluxes 
at  each  detector  position  were  convoluted  with  a RADFLO-calculated 
time  and  wavelength  source  for  a wlde-band  silicon  detector  using  the 
POLY  procedure.  The  detec 'lor  response  function  uoed  in  the  POLY  cal- 
culations is  shown  in  Fig.  9. 

Figure  23  shows  the  source-receiver  geometry  used  in  the  POLO 
procedure.  The  satellite  look  angle,  0 , is  defined  as  the  angle 
between  the  radial  through  the  source  position  and  a line  connecting 
the  subsource  point  and  the  satellite  position. 

The  time- integrated  scattered  fluxes  at  the  satellite  for  a 
ground  albedo  of  0.8  and  look  angles  6 of  0,  45,  75,  and  90  degrees 

la 

are  shown  in  Figs.  24  - 26  and  27  - 29  as  a function  of  source  height 
for  an  atmosphere  containing  a stratus  cloud  layer  and  an  atmosphere 
containing  an  altostratus  cloud  layer,  respectively. 

For  source  altitudes  both  below  and  above  the  cloud,  the  POLO 
procedure  followed  photons  as  they  underwent  collisions  within  the 
cloud.  For  source  altitudes  above  the  clouds,  the  results  of  the 
POLO  calculations  showed  that  the  scattered  flux  at  the  satellite  was 
not  dependent  on  the  magnitude  of  the  ground  albedo.  When  the  source 
was  below  the  cloud,  the  scattered  flux  at  the  satellite  was  highly 
dependent  on  the  magnitude  of  the  ground  albedo. 

The  cumulated  scattered  flux  at  the  satellite  position  as  a 
function  of  retarded  time  for  a source  height  of  0.1  km,  a ground 
albedo  of  0.0  and  a wavelength  of  0.4278  microns  is  shown  in  Fig.  30 
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Fig.  23.  Source-Receiver  Geometry  Uaed  In  Cloud  Problems 
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Fig.  28.  Variation  of  the  Scattered  Flux  at  the  Satellite  with  the 
Height  of  the  Source  for  Look  Angles  of  0,  45,  75,  and  90 
Degrees:  Atmosphere  Containing  Altostratus  Cloud  Layer, 

Ground  Albedo  - 0.8.  X - 0.75  u 
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Fig.  29.  Variation  of  the  Scattered  Flux  at  the  Satellite  with  the 
Height  of  the  Source  for  Look  Angles  of  0,  45,  75,  and  90 
Degrees:  Atmosphere  Containing  Altoatratus  Cloud  Layer, 

Ground  Albedo  “ 0.8,  X - 1.07  u 
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Cumulative  Scattered  Flux  at  Satellite  vs  Time,  Source 


for  each  satellite  look,  angle  when  the  atmosphere  contained  the  low- 
altitude  stratus  cloud.  Retarded  time  Is  defined  as  the  time  after  the 
arrival  of  the  direct  radiation  at  which  the  scattered  radiation  reaches 
the  receiver.  It  Is  seen  that  the  cumulated  scattered  flux  for  a look 
angle  of  0“  Increases  rapidly  with  retarded  time  until  a retarded  time 
of  approximately  10  ysec  after  which  there  Is  very  little  Increase  In 
the  cumulated  flux.  With  an  Increase  In  the  satellite  look  angle,  the 
scattered  flux  starts  to  arrive  at  later  retarded  times.  Similar  results 
were  obtained  for  other  wavelengths  and  for  sources  positioned  below  the 
bottom  of  the  cloud. 

When  the  source  altitude  Is  above  the  top  of  the  cloud  layer, 
then  the  cumulated  distributions  of  the  scattered  flux  vs  retarded  time 
for  look  angles  of  0,  45,  75,  and  90  degrees  and  all  source  wavelengths 
are  similar  to  that  shown  in  Fig.  31  for  the  0.4278-micron  wavelength 
source  at  10-km  altitude  in  the  atmosphere  containing  the  low-altitude 
stratus  cloud  layer. 

The  detectable  power  at  the  satellite  vs  retarded  time  for  times 
up  to  approximately  first  maximum  is  shown  in  Fig.  32  for  a 2-KT  event 
at  0.1-km  altitude  In  the  model  atmosphere  containing  the  stratus  cloud 
layer.  The  curves  shown  are  for  satellite  view  angles  of  0“,  45®,  75°, 
and  90°  and  for  a ground  albedo  of  0.0. 

Increasing  the  source  height  results  in  an  Increase  in  the 
detectable  power  with  the  increase  being  dependent  on  the  satellite  look 
angle.  This  dependence  is  shown  in  Figs.  33,  34,  and  35  for  source 
heights  of  0.5,  3,  and  10  km,  respectively.  Similar  data  for  the  model 
atmosphere  containing  the  altostratus  cloud  layer  are  shown  In  the 
Figs.  36,  37,  38,  and  39  when  the  ground  albedo  is  0.0.  Note  that  the 
curves  for  look  angles  of  0°,  45°,  and  75°  are  much  closer  together 
when  the  source  is  above  the  top  of  the  cloud. 

The  power  that  would  be  measured  at  the  satellite  was  found  to 
be  very  dependent  on  the  magnitude  of  the  ground  albedo  when  the  source 
was  below  the  bottom  of  the  cloud.  When  the  source  was  above  the  top 
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Cumulative  Scattered  Flux  at  Satellite  vs  Time,  Source  Height  “ 10  km 
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of  the  cloud  layer,  it  was  found  that  the  detectable  power  was  not 
dependent  on  the  magnitude  of  the  ground  albedo.  Figure  40  shows  the 
variation  of  the  detectable  power  with  the  ground  albedo  when  the 
source  was  at  0.1-km  altitude  in  the  model  atmosphere  containing  the 
altostratus  cloud  layer.  The  data  shown  are  for  satellite  look  angles 
of  45°  and  90°.  For  the  45°  look  angle,  first  maximum  occurs  at  approxi- 
mately 136  microseconds  in  retarded  time.  Justly  slightly  less  than 
1/2  of  the  power  at  first  maximum  is  reached  at  a retarded  time  of 
70  microseconds  when  the  ground  albedo  is  0.0.  Similarly,  approxi- 
mately 0.1  of  the  power  at  first  maximum  is  reached  at  a retarded  time 
of  36.4  microseconds.  It  is  seen  for  a look  angle  of  45°  that  the 
fraction  of  the  power  detected  at  either  a retarded  time  of  70  micro- 
seconds or  36.4  microseconds  does  not  vary  significantly  with  an  Increase 
in  the  ground  albedo  from  0 to  1.0.  When  the  look  angle  is  increased 
to  90°,  it  is  seen  that  the  variation  of  the  power  at  a retarded  time 
of  192  microseconds  (first  maximum  for  a ground  albedo  of  0.0)  with 
ground  albedo  is  different  from  that  observed  at  retarded  times  of 
106  microseconds  and  76  microseconds  (''^1/2  max  power  and  M).l  max  power). 
The  retarded  time  at  which  first  maximum  occurs  was  found  in  general  to 
Increase  with  an  increase  in  the  magnitude  of  the  ground  albedo.  This 
change  in  the  retarded  time  of  first  maximum  is  also  dependent  on  the 
source  height  and  the  satellite  look  angle. 

Table  XXVI  shows  for  each  cloud  type  the  detected  power  and 
the  retarded  time  at  first  maximum  and  the  retarded  time  at  which  0.1 
the  power  at  first  maximum  is  reached  for  each  source  height  and 
satellite  look  angle.  The  data  shown  are  for  a ground  albedo  of  0.0. 

When  the  source  is  positioned  under  a cloud  layer,  the  retarded  time 
of  first  maximum  tends  to  increase  with  look  angle.  The  retarded  time 
at  which  0.1  the  power  at  first  maximum  is  reached  also  tends  to  Increase 
with  the  look  angle.  When  the  source  is  above  either  cloud  layer  at 
10-km  altitude,  the  retarded  time  at  which  first  maximum  occurs  tends 
to  decrease  with  an  Increase  in  the  look  angle  and  the  time  at  which 
0.1  of  the  power  at  first  maximum  occurs  remains  nearly  constant  for 
all  look  angles. 
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TABLE  XXVI.  POWER  AT  FIRST  MAXIMUM  AND  RETARDED  TIMES  FOR  FIRST  MAXIMUM  AND  WHEN 
THE  POWER  IS  ONE-TENTH  THE  POWER  AT  FIRST  MAXIMUM 


An  attempt  was  made  to  evaluate  the  angular  and  time  distri- 
butions of  the  scattered  light  as  it  exits  from  a cloud  located  above 
a point  isotropic  source.  Computer  runs  were  made  using  the  TPOLOL 
procedure  which  gives  the  time  and  angular  dependence  of  the  scattered 
light  intensity  at  a set  of  point  receivers.  These  runs  were  made  for 
the  10-km  meteorological  range  model  atmosphere  containing  a five  mean- 
free-path  thick  stratus  cloud  located  between  500  and  1000  meters  above 
the  ground  surface.  The  point  source  was  located  at  the  lower  boundary 
of  the  cloud.  An  analysis  of  these  TPOLOL  runs  showed  that  a very 
large  statistical  variation  exists  in  the  results. 

An  example  of  the  statistical  variation  is  illustrated  by 
observing  the  angular  distribution  of  the  scattered  intensity  at  the 
receiver  located  at  a look  angle  of  zero  degrees  (detector  located  at 
the  top  of  the  cloud  directly  above  the  source  position) . Since  the 
source  emission  and  the  atmosphere  are  azimuthally  symmetric  about  the 
source  receiver  axis,  the  scattered  radiation  signal  at  the  detector 
should  also  be  azimuthally  symmetric.  Any  variation  in  the  calculated 
scattered  radiation  signal  with  the  azimuthal  angle  is  then  due  to 
statistical  variations  in  the  results.  Figure  41  shows  the  scattered 
intensity,  integrated  over  time  and  polar  angle,  as  a function  of  azi- 
muthal angle  for  the  receiver  located  at  the  zero-degree  look  angle. 

Although  25,000  histories  were  run  with  an  average  of  37.769  collisions 
per  history,  there  are  variations  by  as  much  as  a factor  of  three  in 
the  azimuthal  distribution. 

It  was  anticipated  that  the  large  statistical  variations  were 
2 

caused  by  the  1/r  factor  in  the  estimating  function.  The  estimate 
of  the  scattered  intensity  contributed  by  each  collision  is  divided  by 
the  square  of  the  distance  from  the  collision  to  the  receiver  position. 

For  a receiver  located  at  the  top  of  the  cloud,  a collision  may  occur 
near  the  top  of  the  cloud  and  thus  produce  extremely  high  estimates  of 
the  scattered  intensity.  A problem,  using  the  model  atmosphere  described 
above,  was  run  with  the  receivers  located  at  a satellite  position 

I 


115 


Detector  Height  = 1 km,  X - 0.4278  y,  5 Mean-Free-Path  Thick  Stratus 


1 


(35800  kilometers  altitude).  The  azimuthal  distribution  of  the  scattered 
intensity  obtained  from  that  problem  (multiplied  by  a factor  of  10  ) is 

plotted  in  Fig.  42.  Note  that  the  variation  in  the  azimuthal  distribu- 
tion for  the  satellite  receiver  is  much  smaller  than  that  computed  for 
the  receiver  at  the  top  of  the  cloud.  This  is  due  to  the  fact  that  all 
collisions  in  the  cloud  and  surrounding  atmosphere  are  approximately  the 

same  relative  distance  from  the  satellite  receiver  and,  therefore,  the 

2 

variation  in  the  1/r  factor  for  each  estimate  of  the  scattered  intensity 
is  small. 

The  statistical  variations  in  the  time  and  spatial  distributions 
of  the  scattered  intensities  computed  for  point  receivers  located  at  the 
cloud  top  were  also  found  to  be  extremely  large.  Therefore,  it  was 
decided  that  area  receivers  rather  than  point  receiver  should  be  used 
to  determine  the  time,  angular,  and  spatial  distributions  of  the  scat- 
tered intensity  at  the  cloud  top.  An  effort  has  been  initiated  to  develop 
a modified  version  of  POLO  that  will  utilize  surface  receivers  Instead  of 
point  receivers  to  compute  the  time,  angular,  and  spatial  distributions 
of  the  scattered  intensity  for  radiation  transport  in  atmosphere  con- 
taining clouds.  It  is  believed  that  the  use  of  surface  receivers  instead 
of  point  receivers  will  reduce  the  statistical  variation  in  the  computed 
scattered  intensities  for  receivers  positioned  within  the  atmosphere. 
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Fig.  42.  Azimuthal  Angle  Dependence  of  the  Scattered  Intensity  at  a Satellite 
Detector  for  a Look  Angle  of  0°:  Ground  Albedo  ■ 0.0,  Source  Height  ' 
500  m,  Detector  Height  “ 35800  km,  X - 0.4278  y,  5 Mean-Free-Path 
Thick  Stratus  Cloud 


VIII.  DISCRETE  ORDINATES  USE  IN  TIME -DEPENDENT 
LIGHT  TRANSPORT  CALCULATIONS 

A study  was  initiated  to  determine  the  feasibility  of  using 
time-dependent  discrete-ordinate  methods  for  light  transport  in  clouds. 

The  extremely  forward  peaked  phase  functions  encountered  in  cloud 
problems  combined  with  the  large  optical  thicknesses  of  clouds  tend  to 
produce  large  statistical  fluctuations  in  results  obtained  with  the 
Monte  Carlo  method.  The  Initial  effort  to  apply  discrete  ordinate 
methods  was  an  attempt  to  utilize  the  TDA  (Time  Dependent  ANISN)  Ref. 

4)  program  distributed  by  the  Radiation  Shielding  Information  Center  at 
OBUIL  as  code  package  CCC-180. 

TDA  is  a one-dimensionax  code,  i.e.,  plane  or  spherical  geometry. 
While  it  is  not  possible  to  describe  realistic  cloudy  atmospheres  con- 
taining a point  isotropic  light  source  with  TDA,  it  is  believed  that  a 
successful  application  of  TDA  to  light  transport  should  precede  applica- 
tions of  two-dimensional  discrete  ordinate  methods  such  as  the  TWO-TIME 
code  currently  undergoing  checkout  at  LASL.  The  attempt  to  utilize  the 
TDA  program  for  light  transport  has  pinpointed  several  difficulties 
that  must  be  overcome  before  discrete-ordinate  methods  can  be  successfully 
applied  to  time-dependent  light  transport  problems. 

In  TDA  it  is  essential  that  the  spatial  intervals  and  time 
boundaries  coincide  such  that  as  the  wave  front  passes  from  one  spatial 
interval  to  the  next,  the  age  (time  after  source  emission)  of  the  wave 
front  is  such  as  to  put  it  in  the  center  of  the  Interval.  Otherwise, 
the  wave  front  becomes  smeared  over  distance  (and  possibly  time)  and 
resolution  is  lost.  For  an  optimum  discrete-ordinates  solution,  the 
spatial  intervals  are  usually  selected  to  be  less  than  one  mean-free-path 
in  width.  After  the  required  spatial  interval  values  are  calculated, 
the  time  required  for  the  wave  front  to  reach  each  spatial  boundary  is 
calculated  using  the  velocity  of  light.  After  the  wave  front  leaves 
the  cloud,  the  spatial  interval  and  time  boundaries  can  be  changed 
(within  reason)  to  suit  the  user's  needs. 
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The  coefficients  of  an  nth  order  Legendre  polynomial  are  required 
to  define  the  scattering  phase  functions  in  the  TDA  program.  The  Legendre 
polynomial  coefficients  for  a given  sized  aerosol  particle  may  be  deter- 
mlnted  with  combinations  of  the  and  scattering  functions  computed 
with  Mie  theory  (Ref.  20)  as  will  be  discussed  later.  However,  the 
expansion  of  the  phase  function  for  aerosols  of  the  sizes  typically 
found  in  clouds  may  require  Legendre  polynomials  of  the  order  of  200  to 
300  or  more.  TDA  requires  that  the  angular  quadrature  set  contain  at 
least  as  many  angle  quadratures  as  there  are  coefficients  in  the  Legendre 
pol)momial  used  to  represent  the  scattering  function.  Use  of  quadratures 
of  this  order  in  the  TDA  program  introduces  problems  in  allocating  core 
space  and  requires  large  amounts  of  computer  time.  It  has  also  been 
determined  that  no  symmetric  Gauss-Legendre  quadrature  would  be  adequate 
(with  a reasonable  number  of  directions)  to  represent  the  peak  portion 
of  the  phase  function. 

Communications  with  one  of  the  authors  of  the  TDA  program  indi- 
cated that  it  may  be  feasible  to  modify  the  TDA  program  to  read  in  a 
precalculated  set  of  quadratures  and  weights.  In  that  case,  a smaller 
quadrature  set  can  possibly  be  selected  to  adequately  represent  the 
scattering  phase  function. 

8.1  Expansion  of  the  MIE  Phase  Function  in  Legendre  Polynomials 

In  applying  results  obtained  from  Mie  theory  to  radiative  transfer 
calculations,  it  is  necessary  when  using  calculational  methods  such  as 
discrete  ordinates  to  expand  the  phase  function,  P(cos0),  in  Legendre 
polynomials : 


P(cos6)  = y a p (cos9)  , 
n“0 

where  p^(cos0)  is  the  Legendre  polynomial  of  order  n. 

Although  it  is  possible  to  calculate  the  coefficients  a^  by 
numerical  integration  from  the  i^^  and  1^  functions  of  Mie  theorj , a 
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direct  calculation  of  these  coefficients  from  theory  is  more  efficient. 
Such  a calculation  scheme  is  provided  in  papers  by  Chu  and  Churchill 
(Ref.  23),  and  by  Clark  et  al.  (Ref.  24).  The  calculations  are  effected 
by  the  use  of  the  Mie  and  scattering  functions,  and  involve 
coefficients  labeled  w and  v which  are  functions  only  of  the  order  n of 
the  Legendre  polynomial  and  are  independent  of  any  other  input. 

Without  discussing  the  methods  of  converting  the  Mie  equations 
to  corresponding  Legendre  polynomial  forms  (see  Refs.  23  and  24),  we  will 
show  and  discuss  the  final  fc’^s  of  the  Legendre  solution  giving  the 
Legendre  coeffic  .ents  a^  by  use  of  computer  programs. 

Equation  (8.1)  gives  the  solution  of  a : 


= Re(AJRe(Aj^)  + lm(A^)Im(A^)  + Re  (B^  )Re  (B^^)  + Im(B^  ) Im  (B^^) 

V = Re(A  )Re(B  ) + Im(A.)Im(B  ) + Re(B.)Re(A  ) + Im(B  )Im(A  ) 
)K  1k  Jk  jk 


(8.1) 

(8.2) 

(8.3) 

(8.4) 


where  Re  = real  part  of,  Im  = imaginary  part  of,  and  A (a,S)  and  B (a, 3) 

n n 

or  (Aj,Aj^,Bj^)  are  Mie  scattering  functions  obtained  with  the  MIE-2 
computer  code  (Rtf.  20). 


The  o::he?  variables  are: 

= 2itR1/X  (8.5) 

RI  = particle  radius  (y) 

A = wave  length  (u) 


3 = a(tn-ik) 

m = index  of  refraction  of  particle  relative  to 
the  surrounding  medium, 


(8.6) 


k = extinction  coefficient  of  the  particle  material. 

*^ST 

K(a,6)=  -^,2  (8.7) 


ST  = scattering  cross  section 
'rr(RI)^  = geometric  geometric  cross  section 


The  definitions  of  w.,  and  v.,  are  given  in  teras  of  the  recursion 

jkn  jkn 

I tionships  used  to  calculate  them  rather  than  the  original  equations 
gi  .n  in  Ref.  23. 


(n+1) 

(8.8) 

V 1.  l,n  ~ 

(2n+lH2n+3)’ 

’ , krt-l,n 

( j +k-  n+1)  ( i +k+n+2 ) 

(8.9) 

(j+k-n+2)  (j+>^-n-*-3)  ’ 

(.i+k+n+1)  (1+k-n)  (;-k+n)  (k-.i+n+] ) 

(8.10) 

w., 

T • 

(j+lc+n)  (j+k-n-1)  (. -k+n-1)  (k-j+n+2)  ’ 

V , ” n(n+l), 

n,  ,n 


(8.11) 


, k-*- 1 , n 
.n 


(j-^k-nr2)  (j+k+n+3 ) 
(.•■-k-n+1)  (j+k+n+2  ’ 


(8.12) 


^j--*,k,n  ^ (j-Mc+n)  (i+k-n-1)  (i+k+n-l)  (k-i+n+2)  (8.13) 

V,.  ^ (j-*-,  -n+l)(i+k-n)(  -k-t-n)(k-j+n+l) 
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A computer  procedure  was  developed  to  evaluate  the  Legendre 
polynomial  coefficients  using  Eq.  (8.1).  The  portion  of  the  MIE-II 
program  (Ref.  20)  which  calculates  the  Mie  scattering  functions  and 
was  incorporated  into  the  program.  The  Legendre  polynomial  coeffi- 
cients computed  for  several  aerosol  sizes  were  used  to  expand  the 
Legendre  polynomials  and  those  expansions  compared  with  the  phase 
functions  computed  for  the  same  sized  aerosols  with  the  MIE-II  program. 
These  comparisons  indicate  the  Legendre  polynomial  coefficients  com- 
puted give  an  accurate  definition  of  the  phase  function.  The  length  of 
the  Legendre  coefficient  array  (a^)  necessary  to  represent  the  phase 
function  appears  to  be  on  the  order  of  2a+2  where  a = 2Trr/X  is  the  aerosol 
size.  Clouds  typically  have  aerosols  up  to  20  to  30  microns  in  radius; 
thus,  for  0.5-mlcron  wavelength  light  the  size  parameter  could  be  as 
large  as  375  requiring  a Legendre  polynomial  of  the  order  of  752. 

Since  the  phase  function  for  an  aerosol  size  distribution 
generally  has  less  structure  than  those  for  the  individual  aerosols, 
the  number  of  combined  coefficients  for  an  aerosol  size  distribution 
possibly  will  not  need  to  be  as  high  as  that  required  for  the  individual 
size  parameters. 
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